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in  using  FLEXSTAB.  Without  their  help,  I doubt  that  this 
paper  would  have  been  completed  in  the  time  allowed. 

Dr.  D.  W.  Breuer,  my  thesis  advisor,  provided  invalua- 
ble assistance  in  analyzing  and  interpreting  much  of  the 
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appreciated. 
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and  proved  to  be  an  excellent  analytical  tool. 
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Abstract 


Using  the  Boeing  FLEXSTAB  digital  computer  system, 
rigid  and  elastic  models  of  the  winglet  configured  KC-135A 
are  made.  With  a rigid  analysis,  the  winglets  reduce  total 
drag  from  2 to  almost  8%  with  improvements  both  laterally 
and  longitudinally  in  static  stability.  Dynamically,  the 
rigid  winglet  model  is  more  stable  laterally  but  slightly 
more  oscillatory  longitudinally.  The  lateral  dutch  roll 
mode  damping  ratio  increased  with  winglets  from  3 to  12% 
with  only  a 3%  increase  in  frequency.  /-The  effect  of  short- 
ening the  vertical  stabilizer  on  the  rigid  winglet  model 
reveals  less  stable  values  for  and  cnr  and  some- 

what less  lateral  dyij1amie-^€ability  than  the  baseline  model. 
For  the  elaa-fic  winglet  model,  the  drag  reductions  noted 
ape'-'fmly  1.3%  due  to  the  "washout"  effect  at  the  wing  tip. 
Elastic  static  stability  is  still  improved  both  laterally 
and  longitudinally  with  winglets;  and  dynamically,  the 
winglets  improve  lateral  stability  with  very  little  effect 
longitudinally.  With  aeroelastic  effects,  the  overall  ben- 
efit derived  from  winglet  application  to  the  model  is  less. 
However,  no  significant  detrimental  effects  due  to  winglets 
are  found. 
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STABILITY  AND  CONTROL  CHARACTERISTICS 


OF  THE  WINGLET  CONFIGURED  KC-135A 


I.  Introduction 


Background 

In  view  of  the  recent  fuel  crisis,  there  has  been 
considerable  renewed  interest  in  ways  to  improve  the 
aerodynamic  efficiency  of  large  aircraft.  Research  in 
this  area  is  critical  when  it  is  realized  that  energy 
use  for  commercial  aviation  alone  has  tripled  in  the  past 
ten  years.  Even  though  fuel  allocations  for  the  most 
part  have  disappeared,  fuel  prices  remain  high,  compris- 
ing, for  example,  as  much  as  40%  of  the  direct  operating 
costs  of  the  commercial  airplane  (Ref  1:19).  Also  there 
is  no  promise  that  fuel  allocations,  and  their  crippling 
effect  on  aviation,  will  not  reappear  in  the  future. 

For  these  reasons,  research  and  development  efforts  are 
currently  being  concentrated  on  ways  to  improve  the  energy 
efficiency  of  existing  and  future  aircraft  designs. 

One  aerodynamic  concept  being  studied  is  the  use  of 
winglets,  which  are  small,  nearly  vertical,  aerodynamic 
surfaces  mounted  at  the  wing  tips.  This  idea  has  been  re- 
viewed extensively  by  Dr.  Richard  T.  Whitcomb  at  the  NASA 
Langley  Research  Center  as  an  effective  method  of  increas- 
ing aerodynamic  efficiency.  Winglets  were  found  to  cause 
a net  reduction  in  induced  drag  by  lessening  the  downwash 
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intensity  of  the  wing  tip  vortices.  NASA  calculations  and 
wind  tunnel  tests  have  shown  potential  for  a 4 to  6%  fuel 
savings  by  current  jet  transports  modified  with  winglets 
{Ref  1:26) . 

A design  study  of  winglets  for  the  KC-135A  and  the 
C-141  aircraft , conducted  by  tha  Boeing  Commercial  Air- 
plane Company  for  the  Air  Force,  showed  a 15.5%  reduction 
in  induced  drag  for  the  KC-135A.  Boeing  also  discovered 
that  winglets  were  more  effective  than  wing  tip  extensions 
by  15.7%  in  reducing  induced  drag  (Ref  2:iii).  This  re- 
sult was  also  verified  by  NASA  Langley  research. 

In  the  NASA  research  program,  a study  was  made  using 
a tailless  full-span  wind  tunnel  model  of  a representa- 
tive jet  transport  aircraft.  Analysis  of  the  model  was 
made  both  with  and  without  winglets.  Their  conclusions 
were  that  winglets  were  stabilizing  both  laterally  and 
longitudinally.  NASA  also  concluded  that  increases  in 
the  directional  stability  of  a winglet-configured  aircraft 
presented  the  possibility  of  reducing  the  vertical  tail 
surface,  thereby  resulting  in  a further  decrease  in  drag 
(Ref  3:7). 

The  U.S.  Air  Force,  in  close  association  with  the 
Boeing  Company,  is  now  giving  serious  attention  to  the 
actual  application  of  winglets  to  an  existing  line  air- 
craft. The  airplane  selected  to  be  retrofitted  and  flight 
tested  for  this  purpose  is  the  Boeing  KC-135A.  The 
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winglet  to  be  used  is  the  one  designed  specifically  for 
this  aircraft  by  the  Boeing  Company.  The  proposed  winglet 
application  is  shown  in  Fig.  1. 

In  view  of  the  approaching  application  phase  of  wing- 
lets  to  the  KC-135A,  considerable  work  must  be  done  in 


teristics  of  the  aircraft  in  its  winglet  configuration. 

The  beneficial  effect  of  induced  drag  reduction  must  be 
compared  carefully  with  any  detrimental  side  effects  that 
may  be  discovered.  Stability  and  control,  aeroelastic 
effects,  and  handling  qualities  predictions  for  both  long- 
itudinal and  lateral  aircraft  motion  are  necessary  before 
the  comprehensive  flight  testing  phase  begins. 

Research  Objectives 

The  primary  objectives  of  this  report  are  to: 

1.  construct  a valid  rigid  and  elastic  mathematical 
model  of  the  Boeing  KC-135A  aircraft  both  with  and  without 
winglets , 

2.  analyze  the  incremental  effect  of  winglets  on 
aircraft  static  and  dynamic  stability  and  control, 

3.  investigate  the  aeroelastic  effects  of  winglets, 

and 

4.  evaluate  the  overall  beneficial  versus  detrimental 
effect  of  winglets. 
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Scope 


Static  and  dynamic  stability  and  control  predictions, 
using  two  rigid  KC-135A  models,  one  with  and  one  without 
winglets,  are  accomplished  for  five  different  flight  con- 
ditions. These  are  specified  in  Table  1.  The  three  cruise 
flight  conditions,  one  heavyweight  climb  condition  and  the 
approach  configuration  with  30  degrees  of  flaps  are  repre- 
sentative of  this  aircraft's  operational  mission  profile 
and  are  chosen  for  this  reason. 

As  an  additional  part  of  the  rigid  stability  and 
control  study,  the  upper  quarter  of  the  aircraft's  verti- 
cal stabilizer  on  the  winglet  model  is  removed  for  flight 
condition  2A.  This  is  done  to  investigate  the  effect  on 
the  lateral-directional  stability  of  the  winglet  configured 
model. 

Aeroelastic  effects  prove  to  be  important  factors  with 
the  addition  of  the  winglet  to  the  wing  tip.  However,, due 
to  time  considerations  and  limited  availability  of  struc- 
tural data  on  the  KC-135A,  the  aeroelastic  analysis  is 
essentially  confined  to  one  basic  aircraft  weight  and  cen- 
ter of  gravity  location.  Nevertheless,  aeroelastic  effects 
on  static  and  dynamic  stability  are  observed  for  three 
different  altitudes  and  dynamic  pressures  as  shown  in 
Table  2. 

With  only  a few  exceptions,  all  of  these  case  stu- 
dies of  the  KC-135A  stability  and  control  are  conducted 


Table  1 

Rigid  Model  Flight  Conditions 


Fliqht  Condition 

Cruise 

1 

Cruise 

2 

Cruise 

2A 

Climb 

3 

Approach 

4 

, 

Altitude  (feet) 

28,500. 

45,600. 

28,500. 

S.L. 

S.L. 

Mach  Number 

.77 

.77 

.77 

.422 

.211 

Gross  Weight 

284,000. 

130,000. 

130,000. 

297,000. 

130,000. 

c.g.  in  % c 

.242 

.321 

.321 

.247 

.321 

Flaps  (degs) 

0 

0 

0 

0 

30° 

y (degs) 

0 

0 

0 

4.01° 

-2.5° 

Table  2 

Elastic  Mr del  Cruise  Parameters 


Flight  Condition 

Cruise  1 

Cruise  2 

Cruise  3 

Altitude  (feet) 

30,000. 

45,000. 

10,000. 

Mach  Number 

.77 

.77 

.422 

q (psf) 

261.67 

128.59 

181.45 

Gross  Weight  268,000  lbs 

Center  of  Gravity  .23  c 
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using  two  models — one  with  and  one  without  winglets.  This 


allows  the  overall  incremental  effect  of  the  winglets  on 


the  aircraft  to  be  observed. 


Analytical  Tools 


The  analytical  tools  that  are  applied  in  satisfying 


the  objectives  of  this  report  consist  of  three  digital  com- 


puter routines — FLEXSTAB,  CALMAT  and  STATVAR. 


FLEXSTAB  is  the  primary  analytical  device  that  is  used 


in  this  report.  It  is  a program  system  which  actually  con- 


sists of  fourteen  separate  programs  of  varying  complexity 
(Ref  4:3-11),  not  all  of  which  are  used  in  this  study. 


FLEXSTAB  has  the  capability  to  evaluate  static  and  dynamic 


stability,  trim  state,  inertial  and  aerodynamic  loading, 


and  elastic  deformations  of  aircraft  flight  configurations 


at  both  supersonic  and  subsonic  speeds.  All  of  these  fea- 


tures of  the  system  are  used  in  this  report;  hov/ever,  only 


for  subsonic  and  near  transonic  flight  regimes.  More  spe- 


cific information  about  FLEXSTAB  is  contained  in  Appendix  D. 


CALMAT  is  an  individual  computer  program  written  for 


this  report  to  calculate  the  aircraft  stability  deriva- 


tives and  control  matrices  for  insertion  into  the  state 


variable  time  histories  program.  CALMAT  also  uses  sub- 


routines to  calculate  the  eigenvalues  and  eigenvectors  of 


the  stability  derivatives  matrix  as  explained  in  Appendix  B. 


STATVAR  is  the  state  variable  time  histories  program 


used  to  analyze  the  aircraft  dynamic  modes  and  responses 
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to  controller  inputs.  Details  concerning  this  program  are 
given  in  Appendix  C. 

Where  needed,  these  digital  computer  programs  are 
executed  using  a Control  Data  Corporation  (CDC)  6600  compu- 
ter with  both  off  line  and  on  line  access  to  a CalComp 
plotter. 

Assumptions 

The  overall  assumptions  made  in  this  report  are  the 
direct  result  of  the  analytical  tools  used  and  the  approach 
to  mathematically  modeling  the  aircraft. 

FLEXSTAB  uses  linear  theories  to  perform  all  of  its 
critical  calculations.  Therefore,  consistent  linear  aero- 
dynamic, structural  and  dynamic  analytical  methods  are  used. 
The  linear  aerodynamic  method  incorporated  is  that  intro- 
duced by  Woodward  (Ref  5) ..  This  method  has  been  extended 
in  FLEXSTAB  to  include  subsonic  flow,  low-frequency  unsteady 
aerodynamics,  and  arbitrary  wing-body-nacelle-tail  arrange- 
ments. Also  viscous  flow  effects  are  neglected  as  with  all 
linear  aerodynamic  potential  flow  programs. 

The  structural  portion  of  FLEXSTAB  used  in  this  study 
consists  of  a structural  elastic  axis  program  that  uses 
beam  theory  to  predict  structural  deformations  due  to  aero- 
dynamic loads  (Ref  4:1).  The  aircraft  fuselage,  and  aero- 
dynamic surfaces  are  divided  up  into  elastic  axis  segments 
which  are  assumed  to  have  constant  stiffness  properties. 

Both  CALMAT  and  STATVAR,  employed  in  the  dynamic  ana- 
lysis , use  the  uncoupled  linearized  aircraft  perturbation 
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equations  of  motion,  with  the  inherent  assumption  that  dis- 
turbances to  the  aircraft  steady  state  flight  condition 
are  small.  This  method  is  commonly  used  in  most  stability 
and  control  problems  and  yields  good  results  when  compared 
to  more  exact  methods.  These  equations  are  developed  and 
discussed  in  detail  in  Appendix  A. 

An  implicit  assumption  in  STATVAR  exists  in  that  the 
program  uses  a numerical  integration  technique  to  predict 
the  aircraft  response  to  control  inputs.  The  assumption  is 
that  the  value  of  the  control  input  over  a sufficiently 
small  time  interval  is  constant. 

Due  to  time  considerations,  much  of  this  analysis  is 
conducted  using  a rigid  airplane  model.  The  rigid  assump- 
tion in  this  case  is  not  very  good;  however,  a constant 
awareness  of  its  limitations  in  analyzing  stability  and 
control  data  does  alleviate  some  of  the  problems  it  causes. 

The  models  usjd  in  this  analysis  have  no  under  wing 

mounted  engine  nacelles  as  does  the  actual  aircraft.  Again, 

this  is  a time  and  also  a computer  cost  consideration  in 

the  modeling  process.  The  incremental  effects  of  the  wing- 

* 

lets  are  of  primary  concern,  and  the  magnitude  of  the  inter- 
related effects  between  the  winglets  and  the  engines  is 
assumed  to  be  negligible.  A ficticious  engine,  however, 
is  modeled  using  the  fuselage  as  a nacelle  in  such  a way 
that  the  thrust  vector,  located  at  the  aircraft  tail,  acts 
along  the  fuselage  reference  line. 


Since  the  principal  aerodynamic  interactions  are  be- 
tween the  wing  and  the  winglet,  thickness  effects  of  the 
vertical  and  horizontal  tail  are  ignored.  The  thickness  of 
the  winglets  and  the  wing  are  modeled  using  airfoil  defini- 
tions from  two  Boeing  documents  .(Ref  2 and  Ref  7). 

In  order  to  evaluate  the  aeroelastic  effects  of  the 
winglets,  the  structural  model  used  is  a modified  model 
taken  from  a FLEXSTAB  case  study  of  the  Boeing  707-320B. 
Again,  this  was  done  primarily  due  to  the  lack  of  pertinent 
structural  information  on  the  KC-135A.  The  707-329B  and 
the  KC-135A  are  geometrically  similar  with  the  exceptions 
of  a longer  fuselage  and  slightly  higher  aspect  ratio  wing 
In  the  case  of  the  707-320B.  The  tail  sections  and  fuse- 
lage contours  of  both  aircraft  are  identical.  Therefore, 
the  assumption  is  that  the  elastic  properties  of  the  two 
aircraft  are  nearly  the  same.  However,  by  shortening  the 
fuselage  and  wing  of  the  707  structural  model  to  make  it 
compatible  with  the  geometry  of  the  KC-135A,  the  result  is 
a stiffer  aircraft.  Also,  on  the  winglet  configured  ela-stic 
model,  the  winglet  is  maintained  rigid. 

Approach  Procedure 

The  first  step  in  beginning  this  evaluation  was  to 
obtain  as  much  geometric,  static  and  dynamic  data  on  the 
Boeing  KC-135A  and  the  Boeing  winglet  as  possible.  This 
information  was  found  in  a Boeing  document  on  -135  Series 
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aircraft  (Ref  7)  and  in  a Boeing  report  on  winglet  design 
and  analysis  (Ref  2) . 

Having  the  necessary  information  on  the  aircraft  and 
the  winglet,  mathematical  models  of  the  KC-135A  with  and 
without  winglets  were  made  using  the  appropriate  programs 
in  FLEXSTAB. 

The  programs  required  in  the  modeling  process  are  the 
Geometry  Definition  (GD) , Aerodynamic  Influence  Coefficient 
(AIC) , Internal  Structural  Influence  Coefficient  (ISIC)  , 
Elastic  Axis  Plot  (EAPLOT) , and  Stability  Derivatives  and 
Static  Stability  (SDSS)  program. 

Using  a scaled  representative  drawing  of  the  KC-135A, 
the  GD  program  is  used  to  calculate  a geometric  description 
of  the  aircraft.  This  data  is  required  for  use  in  almost 
all  of  the  subsequently  used  FLEXSTAB  programs.  The  AIC 
program  uses  the  data  from  the  GD  model  to  calculate  the 
theoretical  aerodynamic  definition.  This  is. done  by  using 
aerodynamic  influence  coefficients  that  satisfy  the  compres- 
sible potential  flow  equations.  The  ISIC  program  calculates 
mass  distributions,  stiffness  and  flexibility  of  the  modi- 
fied 707-320B  structural  model.  Having  the  inertial, 
structural,  aerodynamic  and  geometric  properties  of  the  model 
defined,  the  SDSS  program  predicts  the  rigid  and  aeroelastic 
stability  and  control  characteristics  needed  to  complete 
much  of  this  analysis. 

After  completion  of  the  basic  modeling  process,  the 
model  without  winglets  is  validated  with  the  actual  aircraft 
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by  comparing  aerodynamic  data  and  dynamics  obtained  from 
flight  tests  (Ref  7) . 

Using  the  rigid  models,  which  allow  more  ease  in  chang- 
ing gross  weights  and  center  of  gravity  locations,  the  five 
flight  conditions  of  Table  1 are  compared  with  and  without 
winglets. 

Aeroelastic  effects  are  compared  separately  using 
structural  models.  Elastic  to  rigid  ratios  for  CLa*  cma> 
Cinq*  cni5s*  Cm^,  Cng,  Cyg,  Cip,  and  Cnr  are  calculated 

for  three  values  of  dynamic  pressure;  128.59,  181.69  and 
261.67  psf  and  are  faired  over  a pressure  range  from  0 to 
300  psf.  Deflections  of  the  wing  due  to  the  winglet  and 
the  resulting  effects  on  stability  and  control  are  observed 
for  these  same  dynamic  pressures. 

The  next  section  contains  a detailed  discussion  of 
the  FLEXSTAB  modeling  process  and  a validation  of  the  model 
against  existing  flight  test  data. 
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II.  Validation  of  FLEXSTAB  Computer  Model 
Introduction 

This  section  attempts  to  validate  the  FLEXSTAB  geometry 
and  structural  models  in  order  to  represent  them  as  reason- 
able mathematical  approximations  of  the  actual  aircraft. 
Comparisons  between  the  FLEXSTAB  model  and  flight  test  data 
are  made  for  the  stability  and  control  derivatives  and  the 
longitudinal  and  lateral  dynamic  modes. 

Airplane  and  Winglet  Geometry  Definition 

From  the  aircraft  dimensions  shown  in  reference  7,  the 
KC-135A  geometry  is  modeled  using  the  Geometry  Definition 
(GD)  program  of  FLEXSTAB.  This  program  produces  the  three- 
dimensional  baseline  model  shown  in  Fig.  2.  The  fuselage 
slender  body  and  thin  body  aerodynamic  surfaces  such  as 
the  wing  and  tail  occupy  the  same  relative  locations  as 
those  defined  by  the  actual  aircraft  (Ref  4:31). 

The  panelling  scheme  that  is  used  on  the  model  in 

Fig.  2.  follows  the  summary  of  guidelines  given  in  the 

FLEXSTAB  user's  manual  (Ref  4:47-48). 

¥ 

The  interference  body,  also  shown  with  the  model, 
is  extended  one-half  of  the  root  chord  in  front  of  the 
wing  leading  edge.  Its  purpose  is  to  account  for  the 
upwash  effect  ahead  of  the  wing  attachment  at  the  fuse- 
lage (Ref  4:39). 


A*.  -l  " ' ■■  — 
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Dihedral  angles  on  the  wing  and  the  horizontal  tail 
are  7°  as  shown  on  the  plot.  What  is  not  shown,  however, 
are  the  fuselage  camber  definition,  the  wing  incidence 
angle  and  the  wing  thickness  definition.  These  are  input 
into  the  GD  program  and  are  included  by  FLEXSTAB  where 
necessary  in  all  program  calculations. 

Airfoil  definitions  are  input  into  the  GD  program  at 
five  spanwise  locations  in  order  to  define  the  wing  thick- 
ness. Tail  thicknesses  are  not  defined  on  the  model. 

The  winglet  is  modeled  from  geometry  and  airfoil  de- 
finitions of  the  Boeing  KC-135  winglet  design  of  reference 
2.  Airfoils  are  defined  at  three  spanwise  locations  from 
root  to  tip.  The  winglet  twist  distribution  shown  with  the 
dimensions  in  Fig.  3 is  also  input.  The  dihedral  angle, 
with  reference  to  the  horizontal  plane,  is  70°.  Fig.  4 
shows  the  winglet  configured  geometry  model  of  the  KC-135 
used  throughout  this  report.  A FLEXSTAB  plot  of  only  the 
winglet  is  shown  in  Fig.  5. 

In  order  to  model  the  30°  of  flaps  needed  for  flight 
condition  4 in  the  rigid  analysis,  the  camber  of  the  air- 
foils at  the  respective  flap  locations  along  the  wing  is 

modified.  This  is  done  as  shown  in  Fig.  6.  Realize  that 

» 

this  differs  from  the  actual  aircraft  flap  geometry,  how- 
ever, in  that  the  result  is  a slotted  flap  arrangement. 

The  actual  aircraft  uses  fowler  flaps.  The  differences  that 
result  because  of  this  are  believed  to  be  negligible. 


Fig,  5.  FLEXSTAB  Winglet  Geometry 


pie  Wing  'irfoil  Showing  Camber  Changes  made  for 
Degrees  of  Flaps 


.J  *-* 


The  aircraft  control  surfaces  are  modeled  in  the  Sta- 
bility Derivatives  and  Static  Stability  (SDSS)  program  of 
FLEXSTAB  by  defining  hingelines  and  designating  represen- 
tative panels.  Fig.  7 shows  the  essential  control  surfaces 
for  the  model.  The  entire  horizontal  stabilizer  defines 
the  elevator  in  order  to  observe  more  realistically  the 
trim  conditions  that  result  for  the  various  flight  condi- 
tions. The  rudder  is  defined  with  the  same  geometric  di- 
mensions as  the  actual  aircraft;  but  like  the  vertical 
tail,  it  has  no  thickness.  The  inboard  aileron,  designated 
by  only  one  panel,  is  also  shown.  As  far  as  aileron  control 
is  concerned,  the  KC-135  uses  spoilers  as  well.  However, 
for  the  steady  state  equilibrium  flight  conditions  studied 
here,  the  inboard  aileron  is  believed  to  be  all  that  is 
necessary.  The  spoilers  cannot  be  modeled  using  the  linear- 
ized potential  flow  equations  of  FLEXSTAB. 
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nacelles.  The  707  structural  model  is  modified  by  making 
changes  in  fuselage,  wing  and  engine  configurations. 

Fuselage  modifications  to  the  707  model  amount  to 
shortening  the  fuselage  structure  by  16.83'.  This  is  done 
just  ahead  of  the  wing  root  and  just  aft  of  the  root  trail- 
ing edge.  The  structural  integrity  defined  along  the  wing 
root  and  at  the  tail  is  not  affected  by  this  change.  Also, 
the  shortening  is  accomplished  by  combining  fuselage  elas- 
tic axis  segments  with  nearly  the  same  elastic  properties. 
When  masses  are  involved,  these  are  lumped  together  and 
positioned  in  the  middle  of  the  shortened  interval.  This 
same  procedure  is  used  on  sections  in  front  and  aft  of  the 
wing.  The  center  of  gravity  location  of  the  model  remains 
nearly  unaffected.  However,  the  result  is  a stxffer  fuse- 
lage than  the  707  model,  which,  for  the  shorter  KC-135,  is 
a reasonable  approximation . 

The  wing  is  longer  spanwise  on  the  707  model  by.  5'  over 
the  KC-135.  The  structural  elastic  axis  along  the  entire 
wing  is  therefore  compressed  5'  to  make  it  conform  in  the 
s.ame  proportion  to  the  -135  wing.  In  addition,  this  is  ne- 
cessary in  order  to  align  the  stronger  wing  structure  with 
the  KC-135  chord  segments  where  engines  are  normally  mounted. 
Even  though  the  -135  geometry  model  has  no  wing-mounted 
engines,  there  is  still  a requirement  to  model  the  engine 
support  structure  along  the  wing. 

The  engines  on  the  707  model  create  a different  problem. 


mjK^wa 


is  still  a need  to  represent  the  engine  structural  weight 
effect  on  the  wing.  This  is  done  first  by  numerically 
calculating  the  combined  center  of  gravity  locations  for 
each  engine  and  strut.  Two  rigid  "dumbbells"  of  appropriate 
length  with  total  engine  and  strut  masses  are  then  placed 
at  the  required  elastic  axis  points  on  the  wing.  The  re- 
sult is  shown  in  Fig.  8.  Realize  that  FLEXSTAB  does  not 
require  that  the  aircraft  structural  definition  fall  entirely 
within  the  boundaries  of  the  geometry  model.  Consequently 
the  engine  center  of  gravity  locations  may  be  positioned 
where  required,  which  for  the  KC-135  model  is  in  front  of 
the  wing  leading  edge. 

The  final  result  of  these  modifications  on  the  707-320B 
structural  model  is  shown  in  Fig.  8.  Figures  9 and  10  show 
the  model's  representative  fuselage  and  wing  stiffness 
properties . 

In  adding  the  winglet  to  the  structural  model,  thq 
elastic  axis  of  the  wing  is  extended,  with  the  same  out- 
board stiffness,  to  the  wing  tip.  Three  elastic  segments 
of  nearly  rigid  stiffness  are  then  added  along  the  winglet, 
ending  at  the  nominal  winglet  center  of  gravity  location 
defined  in  reference  2.  Here  the  190  pound  mass  of  the 
winglet  is  added.  Note  that  the  exact  elastic  axis  of  the 
winglet  is  not  modeled  here.  As  a result,  the  axis  used 
is  given  rigid  properties.  Nonetheless,  the  winglet  weight 
effect  on  the  elastic  characteristics  of  the  wing  is 
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Fig.  8.  Baseline  Structural  Model 
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Fuselage  Modex  Stiffness  Distribution 


achieved.  Fig.  11  shows  the  structural  definition  with 
winglets  added. 

For  details  on  FLEXSTAB  internal  structural  modeling, 
see  reference  2. 

Model  Comparisons  With  Flight  Test  Data 

Rigid.  Having  a physical  model  of  the  KC-135  aircraft 
without  winglets,  its  rigid  flight  characteristics  are 
compared  in  Table  3 to  flight  test  data  from  reference  7. 
Note  that  this  flight  test  data  is  semi-elastic  since  the 
longitudinal  derivatives  include  elastic  effects  while  the 
lateral  derivatives  do  not. 

In  studying  Table  3,  realize  also  that  Mach  number, 
altitude,  flight  path  angle,  center  of  gravity  location, 
gross  weight  and  mass  moments  of  inertia  for  the  rigid 
model  are  not  calculated  by  FLEXSTAB  but  are  input  by  the 
user  in  defining  a given  flight  condition.  All  of  the  other 
model  values  shown  are  calculated  by  the  program. 

Trim  conditions,  even  with  the  model's  fictitious  fuse- 
lage reference  line  engine,  compare  fairly  well.  Trim 
angles  of  attack  and  thrust  required  show  the  least  dif- 
ferences in  flight  condition  1.  Nevertheless,  the  other 
flight  conditions  also  show  good  correlation  with  the  flight 
test  trim  values. 

Looking  at  the  longitudinal  data  for  the  f:  !.-;ht  condi- 
tions shown,  only  significant  stability  and  contrw.  deriva- 
tives are  compared.  The  best  correlation  here  in  seen  in  Cl 
cma  and  Cjnq,  however,  are  also  close  in  magnitude. 
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Table  3 

Validation  of  FLEXSTAB  Rigid  Computer  Model 
With  Flight  Test  Data 


Flight  Condition 

3 

4 

(rigid  model) 

^cruise) 

(climb) 

(flaps 

30°) 

Common  Flight  Parameters 

Mach  Number  ! 

.77 

.422 

.211 

Altitude, 

(ft) 

28,500 

S.L. 

S.L. 

Y (degs) 

0 

4.01 

-2.5 

c.g.c 

(sl-ft2) 

.242 

.247 

.321 

I xx  x 106 

2.93 

3.17 

2.05 

Iyy  X 106 

Izz  x 106 

(sl-ft2) 

(sl-ft2) 

4.66  1 

7.48 

4.72 

7.7 

2.46 

4.36 

Ixz  x 106 

(sl-ft2) 

0 

0 

0 

Aircraft  Weight 

284,000 

297,000 

130 

,000 

Longitudinal  Data 


K, 

1 

*F.T. 

FLEX 

*F.T . 

FLEX 

*F.T. 

FLEX 

1, 

6 

®Btrim 

2.40 

2.52 

7.00 

7.07 

-.10 

.07 

5 

P 

Tl  (lbs) 

16,352 

16,323 

38,800 

45,910 

9,000 

13,080 

I- 

CLa  (1/deg) 

.0825 

.0986 

.0715 

.0837 

.0781 

.0791 

1 

!'  1 

Cma  \ 1/deg) 

-.015 

-.0212 

-.012 

-.0173 

-.0175 

-.0106 

I * 

Cinq  (1/rad) 

-14.65 

-18.27 

-13.35 

-15.636 

-14.05 

-x4 .252 

| 

CLl 

.426 

.416 

.463 

.460 

.811 

.804 

1 

CD1 

.024 

.024 

.025 

.039 

.091 

.117 

1 

$ 

Croft  (1/rad) 

-6.57 

-10.22 

-5.47 

-6.79 

-5.52 

-6.23 

1 

♦Flight  Test 

Values 

(Ref  7:6 

.52) 

It 

Lateral 

Data 

! 

Cyg  (1/deg) 

-.0133 

-.0116 

-.0122 

-.0105 

-.0134 

-.0102 

I 

Cs,g  (1/deg) 

-.0035 

-.0025 

-.0031 

-.0022 

-.0040 

-.0022 

1 

h 

Cng  (1/deg) 

.0028 

.0034 

.0024 

.0029 

.0023 

.0026 

k 

§ 

C£p  (1/rad) 

-.3450 

-.4969 

-.3150 

-.4428 

-.3850 

-.4277 

§*  * 

& i 

Cnr  (1/rad) 

-.1940 

-.2926 

-.1720 

-.2648 

-.1860 

-.2580 

i 

CyfiR  (1/deg) 

.0046 

.0076 

.0041 

.0064 

.0040 

.0061 

1 i 

CVdeg) 

.0006 

.0007 

.0005 

.0006 

.0005 

.0005 

1 . 

I 

CnfiH  (1/deg) 

-.0002 

-.0036 

-.0002 

-.0030 

-.0002 

-.0028 

i 

I • 

♦(From  Ref 7: 7 

.4) 
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Laterally  in  Table  3 more  flight  test  derivatives  are 

& i 

s ■ 

available,  and  since  they  are  also  rigid,  closer  compari- 

I 

1 

sons  with  the  FLEXSTAB  values  are  expected.  Note  that  in 

1 

1 

jp» 

general  all  the  flight  test  values  shown  are  reasonably 

li 

S 

close  in  magnitude  to  the  FLEXSTAB  derivatives  with  the 

| 

& 

& 

g 

exception  of  Cn$R* 

s . 

iX 

Table  4 shows  the  longitudinal  and  lateral  dynamic 

i 

! 

% 

period  modes  are  fairly  compatible  with  the  flight  test 

If 

r 

data  shown.  Only  flight  condition  3 shows  a significant 

% 

1 

discrepancy  where  FLEXSTAB  calculates  a phugoid  damping 

i 

ratio  that  is  less  than  half  of  that  for  the  actual  air- 

1  l 
1 ' 

plane.  This  dynamic  mode,  however,  does  not  play  a signi- 

X  ' 

£ . 
& 

s 

ficant  role  in  this  report - 

% 

fi 

& 

£ 

Laterally,  the  model  dynamics  e.o  not  compare  as  close- 

& 

1 

1 

ly.  Dutch  roll  damping  is  consistently  higher  for  the 

1 

model,  and  as  a result,  the  times  to  half  amplitude  shown 

£ 

£ 

1 

are  lower.  Also,  for  all  flight  conditions,  the  model  has 

1 

an  unstable  spiral  mode  whereas  the  flight  test  data  indi- 

1  1 

cates  a stable  spiral  throughout.  The  rolling  mode  shows 

1 I 

It 

the  best  correlation  as  Table  4 indicates. 

1 • 

Elastic.  Table  5 shows  the  model's  mass  moments  of 

1 

! 

inertia  and  center  of  gravity  locations  compared  with  both 

as 

& < 

! * 

the  KC-135A  aircraft  values  and  those  of  the  707-320B 

1 

k: 

structural  model.  Iz2  compares  very  well  with  the  same 

1 

value  for  the  actual  KC-135.  Ixx  and  lyy  do  not  correlate 

f 

as  well,  however.  * 

it 

i i 
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Table  4 

Dynamic  Modes  Comparison 


Flight  Condition  1 

(rigid  model)  (cruise) 


3 

(climb) 


(flaps  30°) 


Short  Period 


wn  (rad/sec) 

period 

(sec) 

^s.p. 

r 

T 1/2 

(sec) 

Phugoid 

wn  (rad/sec) 
period  (sec) 
Sp 

T 1/2  (sec) 
Dutch  Roll 


wn  (rad/sec) 
period  (sec) 
*D.R. 

T 1/2  (sec) 
Spiral 


T 1/2  (sec) 


T 1/2  (sec) 


♦(From  Ref  7:6.55  & 


*F.T. 

FLEX 

■HHI 

FLEX 

*F.T. 

FLEX 

1.67 

1.99 

1.54 

1.81 

1.39 

1.17 

4.07 

3.49 

4.93 

4.17 

6.12 

9.24 

.386 

.418 

.558 

.556 

.676 

.814 

1.07 

.835 

.805 

.687 

.733 

.726 

.0571 

.065 

. .086 

.094 

.157 

.151 

109.7 

96.59 

72.8 

67.13 

39.5 

41.66 

.0339 

.0240 

.028 

.011 

.055 

.080 

355.0 

439.2 

285.0 

663.96 

79.0 

57.08 

% 

1.43 

1.57 

1.31 

1.39 

.95 

.92 

4.4 

4.03 

4.81 

4.59 

6.6 

6.98 

.056 

.129 

.063 

.184 

.069 

.198 

8.6 

3.41 

8.4 

2.71 

10.5 

3.82 

84.6 

unstable 

60.8 

unstable 

44.8 

unstable 

.700 

.558 

.510 

.434 

.520 

.545 

=C.W~- v 4-  - 


■ 


Table  5- 

Mass  Moments  of  Inertia  Comparison 


G.W.  (lbs) 


c.g.  c 

Ixx  x i06  (Si_ft2) 


Iw  x 10°  (sl-ft2) 


1X2  * 106  (Sl-ft2) 
*(From  Ref  7:1. 4-1. 5) 


*KC-135A 

Aircraft 

KC-135A 

Model 

268,500 

268,500 

.25c 

.23c 

2.90 

3.694 

4.60 

3.716 

7.35 

7.256 

- 

.252 

Boeing 

707-320B 


268,500 

.25c 

4.928 

4.966 

9.746 

.245 


Table  6 shows  the  elastic  model  comparison  with  the 
actual  aircraft.  This  is  done  for  elastic  cruise  condition 
1 only.  No  significant  improvement  is  seen  here  over  the 
rigid  model  comparison,  again  shows  the  greatest 

discrepancy. 

Dynamic  modes  for  the  elastic  model  are  not  compared 
because  the  stability  derivatives  show  the  same  relative 
correlation  as  the  rigid  model. 


Conclusion 

Since  the  primary  objective  of  this  report  is  to  inves- 
tigate the  incremental  effect  of  the  winglet  on  the  Boeing 
KC-135A,  the  correlation  shown  in  this  section  between  the 
FLEXSTAB  model  and  the  actual  aircraft  is  adequate. 


The  detailed  analysis  of  the  winglet  effect  on  the 
model  in  the  following  sections  should  do  well  in  predict- 
ing the  trends  in  aircraft  stability  and  control  that  will 
come  out  of  future  flight  tests. 


Table  6 

Elastic  Model  Comparison  With 
Actual  Aircraft 


Cruise  Condition  1 

Common  Flight 

Parameters 

Mach  Number 

.77 

Altitude  (feet) 

30 

,000 

Aircraft  Weight 

268 

,500 

Longitudinal  Data 

F.T. 

FLEXSTAB 

0Btrim  (de9> 

2.4 

2.6 

Ti  (lbs) 

19,536.39 

19,502.00 

CLa  (1/deg) 

.0825 

.0964 

Cma  (1/deg) 

-.0150 

-.0201 

Cmq  (1/rad) 

-14.65 

-17.1325 

cLl 

.421 

.421 

°D1 

.031 

.031 

Cmi  (1/xad) 

-6.57 

-9.1353 

Lateral 

F.T. 

FLEXSTAB 

Cyg  (1/deg) 

-.0133 

-.0137 

C£g  (1/deg) 

-.0034 

-.0024 

Cng  (1/deg) 

.0024 

.0044 

'Ip  (1/rad) 

-.3600 

-.4275 

Cnr  (1/rad) 

-.1940 

-.3524 

cY<5r  (!/ deg) 

.0052 

.0052 

C£fiR  (1/deg) 

.0005 

.0004 

Cn$R  (1/deg) 

-.0002 

-.0025 
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III . Rigid  Model  Analysis  of  the 
Incremental  Effect  o_f  Wingleus 

Introduction 

The  rigid  analysis  to  follow  contains  a detailed  eval- 
uation of  the  winglet  effect  on  drag  and  on  static  and 
dynamic  stability  and  control.  A brief  study  of  the  effect 
of  reducing  the  vertical  tail  on  the  winglet  model  is  also 
made. 

Because  aeroelastic  effects  are  significant,  decisive 
conclusions  from  the  results  of  this  section  alone  cannot 
be  made.  A rigid  study  is  useful,  however,  in  that  it  re- 
veals trends  in  stability  and  control,  trim  states,  and  dy- 
namic behavior  with  the  addition  of  winglets. 

Winglet  Effect  on  Drag  Coefficient  and  Lift  Over  Drag,  L/D 

The  primary  objectives  of  the  Boeing  winglet  design  are 
to  develop  drag  coefficient  reductions  and  subsequent  im- 
provements i".  lift  over  drag  ratios.  Referring  to  Table  7, 
these  objectives  are  satisfied  for  each  of  the  flight  condi- 
tions tested. 

Notice  that  flight  condition  2A  showed  the  greatest 
improvement  in  L/D  and  largest  total  per  cent  decrease  in 
drag.  The  5.58%  increase  noted  in  L/D,  and  the  5.02%  reduc- 
tion in  drag  for  flight  condition  1 compare  well  with  the 
Boeing  wind  tunnel  results  shown  in  Table  8. 
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Table  7 

Winglet  Effect  on  Drag  and  L/D 


Flight  Condition 
(rigid  airplane) 


Mach  Number  M 


Altitude  (feet) 
h 


3 (Flaps  30°) 


% Total  Drag 
Reduction 


.77 

.77  ■ 

28,500 

45,600 

.416 

.427 

5.02 

5.26 

5.58 

5.61 

.422 


.211 


.191  .460  .80- 


7.89  4.09  2.32 


4.34  2.47 


Mach  Number 
M 


Lift  Coefficient 
Cli 


% Total  Drag 
Reduction 


% Increase  in 
Lift/D. a q,  L/D 


Table  8 

Boeing  Wind  Tunnel  Model  Results 
and  FLEXSTAB  Flight  Condition 
1 


Boeing 


.77 


30,000 


.426 


.77 


28,500 


.416 


5.02 


8.4 


5.58 


Static  Steady  State  Stability 

Definitions . As  used  in  this  report,  static  stability 
is  defined  as  the  tendency  of  an  aircraft  to  develop  moments 
or  forces  which  oppose  an  instantaneous  perturbation  of  a 
motion  variable  from  a steady-state  flight  condition.  Steady 
state  motion  is  defined  as  motion  for  which  speed  and  ro- 
tational velocity  remain  constant  with  time  in  a body-fixed 
coordinate  system  (Ref  6:5.1). 

Static  Stability  Criteria.  The  winglet  effect  on  static 
stability  will  be  evaluated  using  the  significant  static 
stability  criteria  shown  in  Table  9.  Also,  mention  of  hand- 
ling qualities  considerations  will  be  made  where  appropriate. 

Longitudinal  Static  Stability.  The  longitudinal  static 
stability  derivatives  discussed  here,  and  shown  in  Table  9, 
are  Ci,a,  and  C^.  In  addition,  the  winglet  effect  on  the 
aircraft  aerodynamic  center  location  and  the  resulting  effect 
on  static  margin  is  studied. 

From  Table  10,  it  can  be  seen  that  a positive  winglet 
effect  on  CLa  is  consistently  noted.  The  implication  is 
that,  with  winglets,  the  aircraft  reaches  a required  lift 
coefficient  at  a smaller  angle  of  attack  than  without  wing- 
lets.  This  is  shown  graphically  in  Figure.  12. 

Of  great  interest  to  static  stability  is  the  position 
of  the  aerodynamic  center  in  its  relation  to  the  aircraft 
center  of  gravity.  As  shown  in  Table  10,  the  winglet 
effect  moved  this  point  aft  in  all  cases  from  7 to  8%. 
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Table  9 

Static  Stability  Criteria  for  Selected 
Longitudinal  and  Lateral  Derivatives 


Perturbation 
Variable  Units 

w (fps) 


a = J1  (degs) 
ui 

6 = —•  (degs) 


Opposing  Force  or 
Moment  Coefficient 


CLa  > 0 


Cma  < 0 


Cnp  > 0 


Table  10 

Winglet  Effect  on  CLa»  xac,  Cma/  and 
Stabilizer  Trim,  6s 


Flight  Condition 
(rigid  airolane) 


% Increase  in  CLot  3.14  3.14  3.14  2.87 


% Aft  Shift  in  xac  8.19  8.3 


7.98  7.23 


(F3ans  30°) 


2.78 


7.08 


% Increase  in 
Static  Margin 

% Negative  In- 
crease  in  Cm« 

% Negative  In- 
crease in  Sstrim 


17.43 

28.13 

27.69 

15.85 

24.03 

21.23 

31.58 

31.01 

19.08 

27.36 

28.5 

38.5 

51.81 

19.62 

46.41 

KG-135A  FLEXSTAB  Model  Flight  Condition  1 


°~No  Winglets 
a ~ With  Winglets 


Angle  of  Attack,  Alpha  (degs) 


Fig.  12.  Winglet  Effect  on  Lift  vs.  Angle  of  Attack 
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Since  the  winglet  increases  tne  effective  angle  of  attack 
of  the  outboard  wing,  this  essentially  moves  the  center  of 
lift  distribution  for  the  total  wing  aft.  This  has  a stabil- 
izing effect  in  increasing  the  aircraft  static  margin  and 
pitch  stiffness,  CmQ: 


Ef; 


i 


i 


i 


1 


i 


i 


Cma  = CLa(xCg  - xac) 

where  (xac  ~ xCg)  i-s  defined  as  the  static  margin.  Table 
10  shows  that  increases  negatively  from  19  to  31%. 

The  stabilizing  changes  in  Cma  and  CLa  noted  thus  far 
with  the  winglet  model  are  significant.  However,  when  viewed 
in  another  way  they  reveal  a slight  trim  penalty.  Referring 
to  Figure  13  for  a .242c  center  of  gravity  location,  the 
winglet  model  shows  a smaller  range  of  useful  trim  lift 
coefficients  and  angles  of  attack  for  a given  range  of  hor- 
izontal stabilizer  settings.  Note  also  that  because  of  the 
more  negative  nose  down  pitching  moment  due  to  angle  of 
attack,  a more  negative  stabilizer  setting  is  required  for 
trim  with  winglets.  This  trend  is  consistent  and  can  also  „ 
be  seen  in  Table  10. 

Lateral  Static  Stability.  The  significant  lateral 
static  stability  derivative  analyzed  here,  and  indicated  in 
Table  9,  is  Cng,  which  defines  the  yaw  stiffness  of  the 
aircraft.  However,  the  winglet  effect  on  Cy^,  and 

Cnr  along  with  Cng  is  also  shown. 

As  indicated  in  Table  11,  the  winglets  cause  a positive 
increment  in  yawing  moment  due  to  sideslip,  Cng.  Note  that 


S 


I 


Angle  of  Attack,  (deg s)  Angle  of  Attack,  (degs) 

Winglet  Effect  on  CLtrim  and  “trim  for  c.g.  = .242c 


Table  11- 

Winglet  Effect  on  Lateral  Stability 
Derivatives 


Flight  Condition 
(rigid  model) 

1 

2 

2A 

3 

4 

(Flans  3C0 

% Negative  Increase 
in  CYB 

14.66 

14.66 

14.66 

14.29 

13.73 

% Positive  Increase 
in  Cnp 

5.88 

6.25 

6.25 

6.90 

3.85 

% Negative  Increase 
in  CiPi 

28.0 

28.0 

30.77 

31.82 

22.73 

% Negative  Increase 
in  C£p 

12.5 

12.38 

12.53 

11.38 

11.06 

% Negative  Increase 
in  Cnr. 

8.03 

7.92 

6.50 

6.57 

5.74 

NOTE:  The  negative  side  force  components  on  the  winglets 

act  behind  and  above  the  aircraft  center  of  gravity. 

a positive  Cng  is  a necessary  condition  in  Military  Airworth- 
iness Requirements  (Ref  6:5.11).  In  comparing  the  incre- 
mental change  in  this  derivative  due  to  winglets,  it  is  small 
but  stabilizing  for  all  five  flight  conditions. 

Observe  in  Table  11  that  all  the  changes  in  the  other 
lateral  derivatives  shown  are  also  stabilizing. 

Summary.  In  the  longitudinal  case,  positive  winglet 
effects  were  observed  in  CLa  and  . The  slight  trim  penalty 
caused  by  Clq  and  was  explained  in  Fig.  13. 

Laterally,  it  was  noted  that  the  winglets  have  a posi- 
tive stabilizing  effect  on  lateral  static  stability. 
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Control  Characteristics  for  Steady  State  Flight 

The  next  area  to  be  analyzed  regards  stability  and 
control  effects  with  winglets.  The  control  aspect  of 
stability  encompasses  two  requirements: 

1.  sufficient  control  power,  and 

2.  no  extraordinary  pilot  effort  (Ref  6:5.19). 
Both  of  these  are  considered  in  the  following  analysis. 

Longitudinal . Here,  the  slight  trim  penalty  with 
winglets  has  already  been  discussed.  As  indicated  in  Table 
10,  the  stabiliser  trim  setting  required  with  winglets  was 
more  negative.  Again,  the  reason  for  this  is  the  more  nose 
down  pitching  moment  effect  of  the  winglets. 

Referring  to  Appendix  E,  Table  I,  the  control  deriva- 
tives Cmfig,  Cl<ss  and  Cd$3  are  the  same  with  or  without  wing- 
lets, which  is  not  surprising.  Apparently,  the  winglets 
are  too  far  outboard  from  the  tail  to  have  any  influence 
on  the  effectiveness  of  the  horizontal  stabilizer. 

Another  practical  consideration  for  longitudinal  stabil 

9 5 g * 

ity  and  control  is  stick  speed  stability,  or  -5^-.  Assuming 
an  unpowered  glide  and  a shallow  glide  angle,  the  uncoupled 
aircraft  longitudinal  equations  become 

in?  * (Cl0  + CLaa  + CL5g5s)qs  and  (1) 

= 0 = Cj^  + cmaa  + cm6s5s 

where  equation  (2)  is  the  moment  trim  condition  (Ref  6:5.21) 
Solving  equations  (1)  and  (2)  for  6S  and  a,  the  result  for 
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r„ ry.-ovr-  .V  -.--t  ^ , . 


6S  can  be  differentiated  with  respect  to  Cl  such  that 


36s  - ~Cmq 

3Cl  CLaCmg  - Cm  Cl* 
s s 


(Ref  6:5.22) 


Since  Cr  = =~,  ~ “ From  this  relation  stick 

qs  3Vp  pSVp3 

speed  stability  becomes  (assuming  constant  Mach  number) 


36s  _ 36s  | 3 Cl  = 4V7  Cma 

"5VP  M ^l|m  pSVp3  Cl“C^  s ^macJL6< 


Since  the  winglet  affect  caused  little  or  no  change  in  Cm*- 

O 

or  CL6s,  the  changes  in  stick  speed  stability  shown  in  Table 
12  are  due  to  changes  in  Cma  and  CLa.  Looking  at  equation 
(4)  in  another  way: 


3 6s  I = 1 

<WplM  cLaCm5s-CL,s< 


where 


(xcjo  ” xac)  S*M* 


Consequently,  equation  (5)  can  be  written  as  a function  of 
static  margin  and  and  Cl^,,: 


36s  1 

"Svp  n “Cmfs  _ 

S.M. 


CL* 

6s 


From  Table  10,  the  winglet  effect  consistently  increases  the 
aircraft  static  margin,  which  results  in  decreasing  the 
denominator  of  equation  (6).  Consequently , the  stick  speed 


Table  12. 

Winglet  Effect  on  Stick  Speed  Stability 


Flight  Condition 
(rigid  airplane) 

% Increase  in 
Stick  Speed 
Stability, 

36e 


1 

2 

2A 

3 

14.29 

23.08 

20.59 

15.48 

(Flaps  30°) 


25.4; 


vH 


stability  increase  of  Table  12,  due  to  the  winglets,  is  veri- 
fied. The  implication  of  this  increase,  from  a pilot's  view- 
point, is  that  more  positive  elevator  deflection  is  required 
to  cause  a change  in  speed  from  a given  trim  speed,  Vpi. 

The  effect  of  stabilizing  increases  in  static  margin  on 
longitudinal  control  is  evident. 

Lateral.  Looking  at  lateral-directional  stability  and 
. vntrol,  the  winglet  had  almost  no  effect  on  the  rudder 
control  derivatives  Cy5R,  Cn5R  and  C£6r,  as  can  be  seen 
in  Table  I of  Appendix  E.  However,  some  small  changes  were 
noted  in  the  inboard  aileron  control  effectiveness  as  shown 
in  Table  13. 

Of  some  concern  here  are  the  FLEXSTAB  calculated  deri- 
vatives of  rolling  moment  due  to  aileron,  Cjt^.  These 
values,  for  both  models  and  all  flight  conditions,  were 
negative.  This  is  not  compatible  with  flight  test  data 
in  reference  7 nor  is  it  reasonable  physically.  The  ex- 
planation for  this  sign  reversal  is  not  understood  and  is 
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rut. 


Table  13 

Winglet  Effect  on  Inboard  Aileron  Lateral 
Control  Effectiveness 


Flight  Condition 

1 

2 

2A 

3 

4 

(rigid  model) 

(Flaps  30°) 

% Decrease  in 

Cy6  (less 
positive) 

29.0 

29.0 

29.0 

25.26 

23.86 

% Decrease  in 

Cl6h  (more 
negative) 

2.75 

2.76 

2.78 

2.08 

1.98 

% Increr.se  in 

Cn6A  (less 
negative) 

3.23 

2.46 

3.78 

2.11 

1.02 

difficult  to  analyze  due  to  the  complexity  of  the  FLEXSTAB 
computer  system.  The  aileron  control  derivatives  Cn,r 
and  cy,$A*  however,  are  acceptable. 

Also  important  to  lateral-directional  stability  and 
control  is  the  capability  to  maintain  straight  and  level 
flight  with  engines  out.  This  is  crucial  to  flying  qual- 
ities requirements  in  that  this  must  be  accomplished  at 
bank  angles  $1  £ 5°  for  all  speeds  larger  than  l-2Vstail 
(Ref  S:para  3.4.12). 

Flight  condition  4 provides  a fairly  realistic  ex- 
ample for  which  engine  out  handling  characteristics  can  be 
evaluated.  In  this  approach  configuration,  the  model  has 
30®  of  flaps  and  is  at  130,000  pounds  gross  weight.  Also 
for  this  condition 


7;f-' 


V = 140  knots  and 


Vstall  = 107  knots 


(From  Ref  9:6-2B) 


. l-2Vstall  = 128.4  knots  < 140  knots 


The  minimum  control  speed,  with  rudder  boost  on,  is  116 
knots  for  a standard  day  at  sea  level  (Ref  10:1A9-14). 
Having  satisfied  the  airspeed  requirements  for  an  engine 
out  approach,  it  is  assumed  that  the  bank  angle  is  also 
less  than  5°.  Using  the  simplified  approach  to  the  lateral 
control  problem  about  the  z axis  only,  the  following  equa- 
tion results  for  zero  sideslip  rudder  deflection: 


6r  = -.rifo- , 
cn<5RqSb 


(Ref  6:5.39) 


Also  assuming  no  pilot  action  at  all 


cngqSb 


(Ref  6:5.39) 


Having  8maxr  an  approximate  way  of  determining  the  amount 
of  lateral  control  needed  to  counteract  any  "sideslip  in- 
duced" rolling  moment  is  available: 


6*  _ "CfcgBmax 
C**A 


(Ref  6:5.39) 


Note  that  equations  (1) , (2)  and  (3)  come  from  the  lateral- 
directional  equations  o...  motion  (Ref  6:5.36).  Using  these 
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simplified  relations,  it  is  possible  to  determine  the 
and  6^  required  v/ith  and  without  winglets  to  counteract  the 


yawing  and  rolling  moments  created  after. engine  failure. 

Instead  of  using  the  questionable  C FLEXSTAB  deri- 
vative, a control  wheel  derivative,  Cj^  , from  reference  7 
is  used.  This  includes,  in  addition  to  the  inboard  aileron, 
the  effect  of  the  outboard  ailerons  and  spoilers,  on  rolling 
moment . 

Using  the  approach  configuration  of  flight  condition 
4 and  the  engine  out  condition  shown  in  Figure  14: 


No  Winglets 


With  Winglets 


Nt  = 150,747  ft-lbs 

C«.  = —.0028  c = .0026 

n«R  nB  -d£5- 


Nt  = 145,814  ft-lbs 

r - -.0028  o _ .0027 
*-nsR  - — f ln»  - — 

**  r?orr  ^ • nPn 


C»B  = -i00|2,*Claw=  ^OOMS 
y deg 


n0  _ —.0027  „ _ .00065 

- ssr>  — asi" 


(q  — 66.04  psf,  S = 2433  ft2,  b = ’0.83  ft.) 


‘(From  Ref  7 :7.4) 


r 

™,l'',|t111’ 111  "“"fc - — > a aj  ~ 

Mb* 

1 

Evaluating  or,  Bmax  and  £w  for  each  model 

| 

No  Winglots 

«■ 

-150,747 

0R  “ (-.0028)  (66.04)  (2433)  (130.83) 

i 

% 

6r  = 2.56°  for  6=0° 

K 

8' 

! 

ir 

g 

-150,747 

8max  “ (.0026)  (66.04)  (2433)  (130.83) 

| 

£ 

I 

b 

Bmax  - -2.76°  , for  fip  = 0 

*? 

S’ 

1 

. _ -(-.0022)  (-2.76) 

w .00065 

6W  - -9.34° 

1! 

With  Winglets 

.-fM  ^ 

-145,814 

°R  ~ (-.0028)  (66.04)  (2433)  (130.83) 

m 

«R  - 2.48°  for  6=0° 

-145,814 

Bmax  = (.0027)  (66.04)  (2433)  (130.83) 

■ /? 

Bmax  = -2.57°  . for  6p  = 0 

ij 

A -(-.0027) (-2.57) 

.00065 

|9 

6„  * -10.68° 

B; 

6r  and  6W  are,  in  both  cases,  within  their  respective  mechan- 

ical  limits.  Note  that  even  though  the  maximum  sideslip 

p 

and  engine  out  yawing  moment  is  less  with  winglets,  the  re- 

1 

guired  control  wheel  deflection  is  slightly  more.  Also, 

1 
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if  ^Rfnax  is  allowed  to  define  the.  minimum  directional 
control  speed,  it  can  be  seen  from  the  above  calculations 
that  this  speed  would  be  slightly  less  for  the  winglet 
model.  This  is  considering  the  same  amount  of  rudder  de- 
flection as  that  required  by  the  baseline  configuration. 

Steady  State  Maneuvering  Flight.  Since  all  the  flight 
conditions  tested  in  this  report  are  for  symmetrical  flight 
only,  the  maneuvering  characteristics  involve  pull  ups  only. 

Based  on  this , FLEXSTAB  calculated  the  elevator  per  g 
derivative  -g— — . This  derivative  is  required  to  be  less  than 
zero  by  handling  quality  requirements  in  reference  8.  It 
has  the  proper  sign  for  both  models  over  the  entire  range 
of  flight  conditions  observed. 

A small  negative  increment  in  elevator  per  g is  consis- 
tently noted  with  the  winglet  model.  Although  this  model 
is  statically  more  stable,  the  implication  here  is  that  it 
is  slightly  less  maneuverable  longitudinally.  This  same 
result  is  also  reflected  in  the  increase  in  stick  fixed 
maneuver  margin  defined  as 


XRef 

3<5e  _ o _ v 

hT  " 0 c9 

> 

4 

See  Table  14  for  changes  in  these  parameters. 

Summary . It  was  seen  that  the  winglets  created  a re- 
quirement for  a negative  increment  in  horizontal  stabilizer 
trim  angles.  The  added  nose  down  pitching  moment  tendency 
from  the  winglets  is  the  primary  cause  of  this  effect  and 
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Table  14 

Winglet  Effect  on  Longitudinal  Control 
Parameters  Steady  Pull-Up 


Flight  Condition 
(rigid  airplane) 

1 

2 

2A 

1 

3 

4 

(Flaos  30°) 

% Aft  Shift  in 
Stick  Fixed  Man- 
euver Point, 

36e  n 

3n 

7.95 

8.03 

7.47 

6.79 

6.25 

% Increase  in 
Maneuver  Margin 

16.00 

24.32 

20.75 

13.30 

14.50 

% Decrease  in 

Longitudinal 

Control/G, 

36e 

3n 

17.41 

25.83 

22.08 

14.47 

15.54 

results  in  a slightly  reduced  Cl  and  a trim  range  for  the 
aircraft.  Also  overall  increases  in  stick  speed  stability 
were  noted.  Laterally,  the  winglets  had  no  effect  on 
rudder  effectiveness  and  only  a small  effect  on  inboard 

* 

aileron  effectiveness.  An  engine  out  example  revealed  less 
yaw  induced  sideslip  and  rudder  required  but  slightly  more 
control  wheel  deflection  with  winglets.  Maneuvering  in  a 
steady  state  pull  up,  the  winglets  were  found  to  offer  in- 
creases in  maneuver  margin  along  with  increases  in  elevator 
required  per  g. 

Dynamic  Stability  Analysis  of  Winglet  Effect 

Introduction.  Dynamic  stability,  as  defined  in  this 
report,  is  the  tendency  of  the  amplitudes  of  the  perturbed 
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motion  of  an  airplane  to  decrease  to  zero  or  to  values 
corresponding  to  a nev;  steady  state,  at  some  time  after 
the  disturbance  has  stopped  (Ref  6:6.1).  Using  this  defi- 
nition, the  effect  of  the  winglets  on  the  dynamic  stability 
of  the  FLEXSTAB  model  is  evaluated.  Note  again  that  the 
dynamic  stability,  studied  in  this  report,  uses  linear  small 
perturbation  theory. 

Because  of  time  considerations,  aircraft  dynamic  modes 
with  and  without  winglets  are  compared  in  detail  for  Flight 
Condition  1 only.  Also  for  this  flight  condition,  sample 
model  responses  to  longitudinal  and  lateral  control  inputs 
are  looked  at.  Here,  changes  in  the  state  variables  are 
plotted  versus  time.  However,  in  an  attempt  to  evaluate 
the  overall  dynamic  effect  of  the  winglets  for  all  of  the 
flight  conditions  tested,  tables  showing  the  winglet  effect 
on  various  dynamic  mode  parameters  are  provided. 

General  Effect  of  Winglets  on  Dynamics . By  way  of  some 

i 

general  remarks,  the  winglets,  when  applied  to  the  rigid 

/ 

model,  provide  no  significant  stabilizing  effer  on  the 
longitudinal  modes.  As  can  be  seen  in  Table  15,  the  de- 
crease in  short  period  damping  ratio  varies  from  over  6 to 
11  per  cent  along  with  an  8 to  almost  13%  increase  in  un- 
damped natural  frequency.  However,  time  to  half  amplitude 
did  decrease  very  slightly. 

The  pitch  damping  derivative,  Qnq , is  important  to 
stability  in  that,  along  with  CLa,  it  determines  short 
period  damping  (Ref  6:5.16).  For  all  the  flight  conditions 
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Table  15. 

Wincjlet  Effect  on  Short  Period  Dynamic 
Mode  and  Cmq 


Plight  Condition 
(rigid  model) 

— T 1 

1 

2 

2a 

3 

3 

(Flaps  30°) 

% Increase  in  Un- 
damped Natural 
Frequency , wn 

— — ... 
9.05 

12.75 

11.81 

8.29 

8.55 

% Decrease  in 
Period,  T 

10.32 

14.87 

17.44 

10.07 

17.69 

% Decrease  in 
Time  to  Half 
Amplitude  T 1/2 

.24 

.44 

.47 

.87 

1.24 

% Decrease  in 
Damping  Ratio 

Gs.P. 

8.37 

11.59 

10.13 

6.65 

6.63 

% Change  of 
Negative  Incre- 
ment in  Cmq 

3.1/ 

2.96 

2.82 

2.85 

2.6 

observed,  the  winglets  resulted  in  a small  stabilizing  change 
in  Cmq,  also  shown  in  Table  15. 

As  far  as  the  phugoid  mode  is  concerned,  almost  no 

f 

change  is  noticed  due  to  the  winglets.  This  long  period 
dynamic  mode  is  influenced  primarily  by  changes  in  CqU/ 
which,  from  Appendix  E,  Table  I,  is  almost  totally  unaffected 
by  adding  winglets  to  the  model. 

Laterally,  the  winglets  are  consistently  stabilizing  to 
the  dutch  roll,  spiral  and  rolling  characteristic  modes. 
Referring  to  Table  16,  a 5 to  almost  14%  decrease  in  time  to 
half  amplitude  and  a 3 to  12%  increase  in  damping  ratio  is 
observed  for  the  dutch  roll  mode.  In  Table  17,  an  average 
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Table  16 

Winglet  Effect  on  Dutch  Roll  Dynamic  Mode 


Flight  Condition 
(rigid  model) 

% Increase  in  Un 
damped  Natural 
Frequency,  wn 

% Decrease  in 
Period,  T 

% Decrease  in 
Time  to  Half 
Amplitude  T 1 

% Increase  in 
Damping  Ratio 

?D.R. 


12.4 


2A 

3 

2.59 

2.88 

2.42 

2.18 

11.44 

9.23 

9.80 

7.07 

(Flaps  30°) 


2. 


Table  17 

Winglet  Effect  on  Rolling  Dynamic  Mode 


Flight  Condition 
(rigid  model) 

% Decrease  in 
Time  to  Half 
Amplitude,  T 1/2 


1 

2 

2A 

3 

9.86 

9.77 

10.43 

10.14 

(Flaps  30°) 


10%  decrease  in  the  time  to  half  amplitude  is  seen  for  the 
rolling  mode.  As  for  the  stabilizing  effect  on  the  spiral 
mode,  the  reader  is  referred  to  Appendix  E,  Table  I. 

Flight  Condition  1^  Dynamics.  Before  analyzing  the  dyna 
mics  of  flight  condition  1,  the  flight  parameters  defining 
it  are  reproduced  below: 
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Altitude  = 28,500  feet 
Mach  = .77 

Gross  Weight  = 284,000  pounds 
c.g.  in  % c = .242 

This  cruise  configuration  has  stability  and  control  deriva- 
tives, from  FLEXSTAB,  tabulated  in  Table  I of  Appendix  E. 

Using  the  FLEXSTAB  data  for  flight  condition  1,  CALMAT 
calculated  longitudinal  and  lateral  stability  derivatives 
matrices,  control  matrices,  and  the  necessary  eigenvalues 
and  eigenvectors,  from  which  STATVAR  time  plots  of  the  cha- 
racteristic dynamic  mode  shapes  were  made. 

Before  going  any  further,  it  is  necessary  to  point  out 
that  the  eigenvalues  calculated  by  CALMAT  are  somewhat  dif- 
ferent than  those  calculated  by  FLEXSTAB.  Table  18  reflects 
these  differences.  The  greatest  variations  between  CALMAT 
and  FLEXSTAB  are  in  the  real  spiral  roots  and  the  real  parts 
of  the  phugoid  roots.  A possible  explanation  for  these  dif- 
ferences is  that  CALMAT  uses  equations  of  motion  that  assume 
cDq/  Cyg,  Cj,£  and  Cng  are  zero.  FLEXSTAB,  on  the  other 
hand,  includes  these  derivatives;  and  as  listed  in  Appendix 
E,  the  computed  FLEXSTAB  values  are  not  ze*>  . Because  these 
modes  have  small  roots  and  are  easily  controlled,  these  diffe- 
rences are  relatively  unimportant.  Note  that  the  short 
period,  dutch  roll  and  rolling  roots  of  CALMAT  and  FLEXSTAB 
compare  quite  well. 

55 


V-*’V  * h.-o ' - <« .> i A\v 1 »' --•»iy V?T~ VA  „ -'  4UA!V'^'  JJ 7^T T 


Table  18 

Comparison  of  CALMAT  and  FLEXSTAB  Eigenvalues 


Flight 
Condition  1 

No  Winglets 

With  Winglets 

Short  Period 

Short  Period 

CALMAT 

(-.8322  + 1. 8051 j ) 

(-.8333  + 2 . 0102 j ) 

FLEXSTAB 

(-.8297  + 1. 8029 j ) 

(-.8321  + 2 . 0080 j ) 

Phugoid 

Phugoid 

CALMAT 

(-.0024  + . 0649 j ) 

(-.0022  + . 0641 j ) 

FLEXSTAB 

(-.0016  + . 0650 j ) 

(-.0015  + . 0642 j ) 

Dutch  Roll 

Dutch  Roll 

CALMAT 

(-.2126  + 1 . 5574 j ) 

(~  2464  + 1 . 6021 j ) 

FLEXSTAB 

(-.2030  + i. 5590 j ) 

(-.2350  + 1 . 6050 j ) 

Spiral 

Spiral 

CALMAT 

(+.0027) 

(+.0009) 

FLEXSTAB 

(+.0044) 

(+.0026) 

Roll 

Roll 

CALMAT 

(-1.2514) 

(-1.3815) 

FLEXSTAB 

(-1.2430) 

(-1.3780) 

As  mentioned  previously,  the  winglet  effect  longitu- 
dinally is  confined  to  the  short  period  mode  of  oscillation, 
with  almost  no  effect  observed  on  the  phugoid  mode.  Figure 
15  shows  the  winglet  effect  on  the  short  period.  It  re- 
veals that  the  winglets  increase  the  relative  amplitudes  of 
the  angle  of  attack  a and  pitch  angle  0 state  variables,  with 
only  small  changes  noted  in  phase  angle  relationships. 

Looking  at  time  history  plots  in  Figures  IGa  and  16b  from 
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STATVAR,  tho  amplitude  differences  in  a and  0 are  also  appa- 
rent. Almost  no  differences,  however,  in  amplitudes  or 
phase  relations  are  seen  in  u and  q.  The  overall  conclusion 
is  that  the  effect  of  the  winglets  on  the  short  period  mode 
is  small.  Table  19  shows  the  short  period  dynamic  charac- 
teristics. 

Laterally,  the  winglet  effect  on  the  dutch  roll  mode 
is  the  most  important.  The  spiral  and  rolling  modes  are 
nonoscillatory  and  are  usually  well  within  the  pilot's  cap- 
ability to  control.  Suffice  it  to  say  that  the  winglets 
are  stabilizing  to  these  modes  as  can  be  seen  in  Table  20. 

The  derivative  Cnr  usually  has  the  most  effect  on  spiral 
stability  while  C£p  primarily  affects  the  rolling  mode  (Ref 
6:6.77-6.85).  The  winglet  effect  on  these  derivatives  is 
also  shown  in  Table  20. 

As  previously  noted,  the  winglet  effect  on  the  dutch 
roll  lateral  aircraft  mode  is  stabilizing.  For  flight  con- 
dition 1,  the  overall  effect  of  the  winglets  on  this  dynamic 
mode  is  shown  in  Table  21.  From  this,  it  can  be  seen  that 
the  directional  stability  derivative  Cng,  which  normally  has 
a strong  effect  on  dutch  roll  frequency  and  damping,  changes 
little  with  winglets.  Consequently,  the  dynamics  shown  in 
Table  21  are  mostly  affected  by  changes  in  Cnr  and  Cy^  (Ref 
8:134).  These  derivatives  have  a powerful  effect  particular- 
ly on  dutch  roll  damping  (Rei  6:6.77-6.85). 

The  dutch  roll  mode  shapes  for  flight  condition  1 are 
shown  in  Fig.  17.  Here  it  can  be  seen  that  little  change 
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Table  19 

Winglet  Effect  on  Short  Period  Dynamic  Behavior 


Flight  Condition  1 
(rigid  model) 

! 

No  Winglets 

With  Winglets 

wns.p.  (rad/sec) 

1.99 

2.17 

Cs.p. 

.418 

.383 

Ts.p. 

3.49 

3.13 

Ts. p.1/2  (sec) 

.835 

.833 

Table  20 

Winglet  Effect  on  Spiral  and  Rolling  Modes 


Flight  Condition  1 
(rigid  models) 

i 

No  Winglets 

With  Winglets 

Spiral  Divergence 
T2  (secs) 

156.35 

263.42 

Rolling  Time  to 
Half  Amplitude 

.558 

.503 

C£p  (1/rad) 

-.4969 

-.5590 

Cnr  (1/rad) 

-.2926 

-.3161 

in  the  state  variable  relationships  results.  The  amplitudes 

of  jJ  and  9 are  reduced  slightly  along  with  small  changes  in 

r and  f . Also  phase  angle  relationships  remain  fairly 

constant  with  or  without  winglets.  Figures  18a  and  18b 

% 

show  the  dutch  roll  state  variable  response  versus  time 
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Table  21 

V7inglet  Effect  on  Dutch  Roll  Dynamic  Behavior  and 
Selected  Lateral  Derivatives 


Flight  Condition  1 
(rigid  model) 

I 

| No  Winglets 

V7ith  V7inqlets 

ttno.R.  (rad/sec) 

1.57 

1.62 

Cd.r. 

.129 

.145 

tD.R.  (sec) 

4.03 

3.91 

Td.R.1/2  (sec) 

3.414 

2.949 

Cnr  (1/rad) 

-.2926 

-.3161 

Cyg  (1/deg) 

-.0116 

-.0133 

Cng  (1/deg) 

| 

.0034 

i 

.0036 

using  the  real  parts  of  the  eigenvectors  shown  in  Fig.  17 
as  initial  conditions.  Here  the  amplitude  differences  are 
more  easily  seen.  The  greatest  apparent  winglet  effect  is 
in  reducing  the  amplitude  of  0. 

Figures  19  (a  & b)  and  20  (a  & b)  attempt  to  show  the 
winglet  effect  on  the  longitudinal  and  lateral  state  varia- 
bles for  two  control  inputs — a 20°  elevator  impulse  and  a 
10°  rudder  impulse. 

In  Fig.  19  (a  & b) , the  model's  response  to  a 20°  ele- 
vator impulse  shows  only  small  differences  in  the  changes 
of  the  state  variables  compared  with  and  without  winglets.* 


♦Note  that,  20°  of  elevator  here  is  actually  20°  of  stab- 
ilizer which  exceeds  the  mechanical  limits  on  the  actual  air- 
craft. The  stabilizer  was  used  for  amplification  purposes. 
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Histpry  Plots  of  Dutch  Roll  With  and  Without  Winglets 


16.00  24.00  32.00  0.00  8.00  16.00  24.00  32.00 

SECONDS)  . T (SECONDS) 

/ 

Winglets  With  Winglets 


The  most  noticeable  difference  is  in  u.  Its  peak  value  is 
slightly  greater  without  winglets.  Also,  to  a lesser  ex- 
tent, the  amplitudes  of  the  other  variables  are  less  with 
winglets  although  this  is  difficult  to  see. 

Laterally,  for  impulsive  10  degree  rudder  inputs,  the 
winglets  do  not  appear  to  affect  the  frequency  or  magnitudes 
of  the  state  variables.  This  is  shown  in  Fig.  20  (a  & b) . 

Summary.  The  overall  winglet  effect  on  the  longitudi- 
nal dynamic  modes  was  felt  primarily  in  the  short  period. 
Here  decreases  in  damping  ratios  varied  from  over  6%  to  al- 
most 12%.  Laterally,  the  dutch  roll  mode  was  more  stable 
with  winglets.  Damping  ratios  increased  from  3 to  12%. 
Positive  effects  on  the  rolling  and  spiral  modes  were  seen 
also.  As  far  as  the  model  responses  to  control  inputs  are 
concerned,  only  very  small  changes  in  the  state  variable 
responses  were  seen  with  winglets  added  to  the  model. 

Analysis  of  the  Effect  of  a Shorter  Vertical  Tail  on  the 
winglet~ Mo9eT 

Introduction.  As  it  affects  lateral-directional  and 
dynamic  stability  and  control,  the  purpose  of  this  section 
is  to  investigate  the  possibility  of  reducing  the  vertical 
stabilizer  on  the  wir.  glet  configured  aircraft  model  as 
shown  in  Fig.  21.  No  reduction  of  the  rudder  control  sur- 
face itself  is  involved. 

Because  the  shortened  \ ertical  stabilizer  has  no  thick- 
ness definition,  FLEXSTAL,  in  its  calculations,  shows  no 
effect  on  the  longitudinal  characte’ Istics  of  the  model. 

'Ml 


Also,  no  observation  of  the  additional  drag  reduction  in- 
volved can  be  made.  In  reality,  such  a reduction  would  be 
small,  and  the  primary  interest  here  is  in  lateral  stability 
and  control.  It  has  been  proposed  that  the  added  lateral 
stability  provided  by  the  winglets  could  make  the  shortened 
vertical  tail  configuration  of  Fig.  21  feasible.  The  com- 
plete data  calculated  by  FLEXSTAB  for  this  model  is  shown 
in  Table  II  of  Appendix  E. 

In  making  this  analysis,  the  lateral  characteristics 
of  the  short  tail  configuration  are  compared  with  those  of 
the  full  tail  model  with  and  without  winglets.  This  is  done 
only  for  flight  condition  2A,  defined  as  cruise  equilibrium 
flight  at  28,500  feet  c.ltitude,  .77  Mach  number  and  130,000 
pounds  gross  weight. 

Lateral -Directional  Stability.  Here  selected  lateral 
stability  derivatives  are  discussed.  These  are  Cyg,  Cng, 
CiQt  C*p  and  Cnr. 

From  Table  22,  it  can  be  seen  that  the  side  force  deri- 
vative, Cyg,  becomes  less  negative  with  the  short  tail 
winglet  model.  This  is  not  surprising  due  to  the  strong 
effect  the  vertical  tail  surface  has  on  this  derivative. 
However,  Cyg  with  this  configuration,  is  a less  stable  value 
than  that  calculated  for  the  full  tail  model  without  wing- 
lets by  4.3%.  Also,  when  compared  to  the  full  tail  winglet 
model,  this  difference  in  Cyg  is  17.u6%. 

The  yaw  stiffness  derivative,  Cng,  moves  in  a less 
stable  direction  by  shortening  the  vertical  tail.  When 
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Table  22 

Lateral  Derivatives  Comparison  With 
Short  Tail  Winglet  Model 


Flight  Condition  2A 
(rigid  models) 

N.W. 

(full  tail) 

W.W. 

(full  tail) 

W.W. 

(short  tail) 

Cyg  (1/deg) 

-.0116 

-.0134 

* -.0111 

Cng  (1/deg) 

.0032 

.0034 

* .0024 

C*6  (1/deg) 

-.0026 

-.0034 

-.0029 

Cip  (1/rad) 

-.4959 

-.5577 

-.5533 

Cnr  (1/rad) 

-.2677 

-.2851 

* -.2237 

*Note  that  these  ’:<aiues  are  less  stable  than  those  for  the 
no  winglet  model  in  the  first  column. 

compared  to  the  full  tail  winglet  model,  Cng  is  29.41%  less 
stable.  Cng  is  a less  stable  value  than  that  for  the  full 
tail  model  without  winglets  by  25%. 

Cj,g.,  or  the  dihedral  effect  derivative,  becomes  14.71% 
less  stable  than  the  same  value  for  the  full  tail  v-inglet 
:;odel  shown  in  Table  22.  Notice  that  here  the  winglets  are 
able  to  create  a compensating  effect  on  this  derivative  to 
make  it  11.54%  more  stable  than  the  value  calculated  for  the 
baseline  modex. 

The  destabilizing  effect  of  reducing  the  vertical  tail 
had  the  least  effect  on  C£p.  Again,  the  winglets  still 
have  a strong  compensating  influence  on  this  derivative  by 
making  it  11.57%  more  stable  than  the  value  show.t  for  the 
no  winglet  model. 
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Cnr/  or  the  yaw  damping  derivative,  is  not  as  stable 
with  the  shorter  tail  when  compared  to  either  of  the  full 
tail  models.  Here  again,  the  winglets  fail  to  offer  a com- 
pensating effect.  Cnr  for  the  short  tail  winglet  model  is 
16.44%  less  stable  than  the  same  value  for  the  baseline 
configuration. 

Lateral  Control.  Because  of  the  12%  reduction  of  the 
vertical  tail  area,  the  result  is  a less  effective  aero- 
dynamic surface.  The  value  of  Sv,  the  planform  vertical 
tail  area,  is  smaller.  This  effectively  reduces  the  rudder 
control  authority  as  Table  23  indicates.  Note  that 


CyfiR  = CL6va6Rnv  ' 


c’  • < - CL«va5K'"'  T IT 


(Ref  6:4.44) 


(Ref  6:4,40) 


-cLavasRnv  If  ^ 


(Ref  6:4.79) 


These  equations  show  the  linear  dependence  of  CygR,  C£$r, 
and  on  Sv.  For  a given  dynamic  pressure,  nv  = qv  ; 

q 

the  rudder  deflection  required  in  controlling  a yawing  mo- 
ment, rolling  moment  and  side  force  is  increased.  Note 
that  the  percentages  shown  in  Table  23  are  independent  of 


Table  23 

Short  Tail  Effect  in  Reducing  Rudder  Authority 
f°r  Cy5R,  CnSR  and  Ch4r 


Flight  Condition  2A 
(rigid  mode 1 s ) 


Compared  to  Full  Tail  Models 
With  or  Without  Winqlets 


*%  Decrease  in  C, 
(less  positive) 

% Decrease  in  C 
(less  positive) 


21.43% 


*%  Increase  in  CT 
(less  negative) 


11.76% 


*Note  that  these  reductions  in  rudder  effectiveness  in  side- 
force  and  yawing  moment  are  comparable  to  the  12%  reduction 
of  the  vertical  tail. 


winglets.  It  was  noted  earlier  that  the  winglets  themselves 
have  almost  no  effect  on  rudder  authority. 

Dynamic  Lateral-Directional  Stability.  Table  24  shows 
the  comparison  of  the  short  tail  model  dynamics  with  both 
of  the  full  tail  model  configurations.  Hers  the  effect  of 
the  shorter  vertical  tail  is  primarily  on  the  dutch  roll  and 
spiral  modes.  No  effect  is  noted  on  the  rolling  mode. 

Looking  at  the  dutch  roll  mode  in  Table  24 , the  short 
tail  model  has  a longer  time  to  half  amplitude  and  a smaller 
undamped  natural  frequency  than  either  of  the  other  two 
full  tail  models.  Damping  ratio,  however,  falls  between 
the  value  for  the  baseline  model  and  the  full  tail  winglet 
model.  This  value  is  3.27%  better  than  the  baseline 
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Table  ?4 

Dynamic  Modes  Comparison  With 
Short  Tail  Winglet  Model 


Flight  Condition  2A 
(rigid  models) 

N . W . 

(full  tail) 

W.W. 

(full  tail) 

W.W. 

(short  tail 

Dutch  Roll 

wn  (rad/sec) 

1.93 

1.98 

1.66 

period  (sec) 

3.3 

3.22 

3.83 

?D.R. 

.153 

.168 

.158 

T 1/2  (sec) 

2.36 

| 

2.09 

| 

2.65 

Spiral  (FLEXSTAB) 

tn 

II 

i 

• 

o 

o 

u> 

(S=~. 0049) 

(S=- . 0064 ) 

T 1/10  (sec) 

675.68 

470.42 

358.3 

T-  1/2  (sec) 

203.4 

141.61 

107.86 

Spiral  Approximation 

(S=-. 0604) 

(S=-. 109) 

(S=-.1044) 

T 1/10  (sec) 

38.12 

21.18 

22.06 

T 1/2  (sec) 

11.48 

6.36 

6.64 

Roll 

T 1/10  (sec) 

1.242 

1.113 

1.113 

T 1/10  (sec) 

.374 

.335 

.335 

configuration.  When  comparing  the  full  tail  and  short  tail 
winglet  models,  there  is  a 5,95%  reduction  in  damping  ratio 
caused  by  the  shorter  vertical  tail.  Tne  overall  net  effect 
is  a less  stable  dutch  roll  mode. 


•B*  jLA!<6fcM~  -■, 
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As  far  as  the  spiral  mode  is.  concerned,  the  FLEXSTAB 
values  in  Table  24  show  a large  increase  in  spiral  stabil- 
ity for  the  short  tail  winglet  model  over  both  of  the  full 
tail  configurations.  This  is  an  unexpe  ted  result,  parti- 
cularly after  having  noted  the  consistently  less  stable 
derivatives  discussed  for  this  model  earlier.  Using  the 
spiral  approximation  shown  in  reference  6,  page  6.54,  a 
more  reasonable  trend  is  found  in  comparing  the  three  models. 
Because  this  approximation  ignores  the  angle  of  bank  degree 
of  freedom,  the  large  differences  in  the  magnitudes  of  the 
time  characteristics  compared  to  the  FLEXSTAB  values  are 
somewhat  expected.  However,  the  important  point  of  the 
comparison  is  that  the  spiral  for  the  short  tail  winglet 
model  is  less  stable  than  the  full  tail  winglet  model,  which 
is  a more  intuitive  result.  Without  a more  detailed  ana- 
lysis, the  contrast  among  the  FLEXSTAB  calculated  spiral 
mode  characteristics  is  not  understood. 

Summary . Reducing  the  vertical  tail  by  12%  on  the  wing- 
let model  results  in  less  stable  values  for  Cyg,  Cng,  Cj£g, 

Cj,p  and  Cnr.  Cyg,  Cng,  and  Cnr  were  noted  to  be  -ass  sta- 
ble than  those  values  calculated  for  the  full  tail  baseline 
model. 

As  for  lateral  control  effectiveness,  the  shorter  tail, 
and  hence  smaller  Sv,  revealed  reductions  in  rudder  authority 
for  Cy$R,  ar.d  CngR.  This  was  verified  by  equations  from 
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Dynamically,  the  dutch  roll  mode  for  the  shortened 
vertical  tail  model  was  less  stable  than  either  of  the  two 
full  tail  models.  Improvements  in  spiral  stability,  with 
the  short  tail  winglet  model,  actually  surpassed  that  of 
the  full  tail  winglet  model.  This  result  was  very  surpris- 
ing and  is  not  understood.  If  one  can  accept  the  • >pari- 
son  made  with  the  spiral  approximation,  the  FLEXSTAB  results 
seem  questionable.  Finally,  the  rolling  mode  remained  un- 
affected by  the  change  in  the  vertical  tail. 


IV.  Elastic  Model  Analysis 

Introduction 

As  mentioned  earlier,  aercelastic  effects,  particularly 
with  large  aircraft  such  as  the  KC-135A,  are  significant. 
Furthermore,  the  weight  of  the  winglet  and  its  aerodynamic 
forces  and  moments  create  additional  elastic  deflections 
of  the  wing.  These  deflections  yield  somewhat  different 
winglet  effects  than  the  rigid  model  analysis  of  the  pre- 
vious section. 

Winglet  Effect  on  Drag  Coefficient  and  Lift  Over  Drag  L/D 
As  with  the  rigid  analysis,  the  winglet  produces  drag 
reductions  and  improvements  in  L/D  ratios.  However,  using 
the  elastic  models,  these  improvements  are  much  less  as  shown 
in  Table  25.  For  all  the  cruise  conditions  observed,  drag 
reductions  due  to  winglets  amount  to  only  1.3%.  Consequently, 
L/D  ratios  are  not  greatly  improved  either.  In  reducing 
drag,  the  elastic  winglet  model  is  4 to  6 per  cent  below 
the  figures  for  the  rigid  model. 

Winglet  Elastic  Effect  at  the  Wing  Tip 

In  order  to  place  the  following  part  of  this  section 
into  a proper  physical  perspective,  some  discussion  of  the 
elastic  winglet  effect  on  the  wing  is  'ded. 

Table  26  shows  the  "washout"  angles  and  fuselage 
reference  line  angles  of  attack  with  and  without  winglets. 

79 


Table  25 

Winglet  Effect  on  Dray  and  Lift/Drag,  L/D 


Cruise  Condition 
(static  elastic  models) 


q (psf) 


Mach  Number,  M 


Altitude  (feet) , h 


Airplane  Lift 
Coefficient,  Cl 


% Increase  in  L/ 


% Decrease  in  Drag 


181.69 


Table  26 

Washout"  Angles  and  Angles  of  Attack 


Cruise  Condition 
(static  elastic  models) 


(q=261.67psf)  I (a=128 . 57psf ) (q=181 . 69psf ) 


-3. 


Washout  Angle  at  the 
Wing  Tip  (degs) 


Angle  of  Attack 
of  Fuselage  (degs) 


-2.9 

-3.9 

-2.9 

1 

U> 

• 

CO 

7.262 

7.212 

5.716 

5.687  1 

"Washout,  as  used  here,  means  a negative  rotation  about 
the  structural  y-axis  of  the  wing,  i.e.,  leacing  edge  down 
Note  that  the  model  without  winglets  displays  a fairly 


constant  -3°  "washout"  at  the  tip.  The  winglet  model  has 
a degree  more  "washout."  This  effect  is  due  to  three  things: 

1.  a small  spanwise  shift  in  the  pressure  distribu- 
tion outboard  along  the  aft  chord  section  of  the 
wing, 

2.  an  aft  chordwise  shift  in  pressure  at  the  out- 
board wing,  and 

3.  the  aerodynamic  loads  created  on  the  winglet. 

All  three  of  these  phenomena  create  forces  which  add  an  in- 
crement in  negative  torque  at  the  wing  tip.  Fig.  22  demon- 
strates, graphically,  the  aft  shift  in  pressure  at  the  wing 
tip  for  cruise  condition  1.  This  effect  is  not  seen  in- 
board on  the  wing.  In  fact,  the  chordwise  pressure  distri- 
butions on  the  inboard  wing  sections  are  the  same  with  or 
without  winglets.  Also  in  Fig.  22  there  is  a small  spanwise 
movement  of  pressure  implied  along  the  aft  chord  section 
of  the  wing.  Average  pressures  for  each  chord  section  were 
taken  and  plotted  along  spanwise  locations  of  the  wing  in 
an  attempt  to  show  this  over  the  whole  wing.-  However,  this 
effect  is  found  to  be  very  small  and  confined  only  to  the 
most  outboard  section  of  wing  panels.  It  should  be  pointed 
out  that  pressure  plots  of  the  outboard  wing  are  somewhat 
degraded  by  neglecting  the  viscous  interaction  at  the 
juncture  of  the  wing  and  the  winglet. 

Perhaps  the  greatest  contribution  to  the  increase  in 
"washout"  at  the  tip  comes  from  the  aerodynamic  forces  on 
the  winglet.  Since  the  vortex  diffusing  effect  of  the 
winglet  results  in  a more  lift  effective  outboard  wing, 
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CRUISE  CONDITION  1 


Fig.  22.  Chordwise  Pressure  Distribution  of  Outboard  Wing 
With  and  Without  Winglets  (elastic  modals) 
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the  elastic  winglet  model  experiences  greater  bending  and 
consequently  more  washout  due  to  the  sweep  of  the  wing. 

The  sweep  effect  is  shown  in  Fig.  23. 

Referring  to  the  designated  chord  sections  or.  the 
FLEXSTAB  winglet  plot  of  Fig.  24,  winglet  chordwise  pres- 
sure distributions  for  a rigid  flight  condition  sample  of 
the  previous  section  are  shown  in  Fig.  25.  As  expected, 
the  pressure  is  generally  higher  at  the  leading  edge  then 
tapers  off  moving  aft  for  each  chord  section  shown.  Looking 
at  Fig.  26  for  the  winglet  on  the  elastic  model,  however, 
the  reverse  is  shown.  There  is  a larger  pressure  rise 
indicated  toward  the  trailing  edge  of  tha  winglet.  This 
effect,  and  the  fact  that  for  an  airfoil  the  pressure  dif- 
ferential is  normally  highest  at  the  leading  edge,  make 
these  results  questionable.  No  explanation  for  these 
elastic  winglet  pressure  distributions  is  available.  The 
winglet  pressure  plots  on  the  rigid  model,  however,  seem 
credible. 

Referring  again  to  Table  26,  observe  the  almost  negli- 
gible winglet  effect  on  angle  of  attack.  This  is  a "washout" 
effect  which  was  absent  in  the  rigid  analysis.  Where  an 
angle  of  attack  reduction  was  consistently  observed  for 
the  rigid  winglet  model  because  of  a more  lift  effective 
outboard  wing,  the  elastic  winglet  model's  increased 
"washout"  tendency  reduces  this  effect.  Nonetheless,  the 
winglet  still  offers  some  vortex  diffusing  effects  at  the 
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Fig*  26.  Winglet  Pressure  Distribution  Plots  for  Elastic 
Model 
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Static  Steady  State  Stabilitv 


Using  the  definitions  and  criteria  for  static  steady 


state  stability  defined  in  the  previous  section,  the  elas- 


tic winglet  model  is  evaluated  longitudinally  and  laterally 


by  comparison  with  the  baseline  model. 


Longitudinal  Static  Stability.  Referring  to  Table  27, 


all  cruise  conditions  show  an  aft  shift  in  aerodynamic 


center  location  due  to  winglets.  These  figures  are  roughly 


3 to  4%  less  than  those  noted  for  the  rigid  model.  Also, 


important  to  point  out  in  Table  27,  is  that  the  static 


margin  increased  7 to  11%.  This  has  a positive  stabilizing 


effect  on  Cma  from  equations  previously  related.  Realize 
that  this  increase  in  static  margin  is  little  affected  by 


the  winglet  induced  aft  shift  in  the  overall  center  of 


gravity  location  of  the  model.  This  value  is  only  -f. 0022c 


aft.  Consequently,  any  weight  penalty  effect  due  to  the 


winglet  is  negligible  in  this  case. 


From  Table  28  there  is  a large  difference  between  the 


rigid  Cl0  increase  with  winglets  and  that  noted  using  the 


elastic  model.  In  fact,  the  elastic  increases  are  almost 


negligible.  This  degradation  in  the  elastic  winglet  effect 


on  CLa  is  again  due  largely  to  the  "washout"  effect  at  the 


Having  noted  increases  in  both  static  margin  and  CLa, 


Cma  becomes  more  negative  by  the  amounts  shown  in  Table  28. 


Table  27 

Winglet  Effect  on  Aerodynamic  Center  Location 


Cruise  Condition 
(static  elastic  models) 


% Aft  shift  in  Aero- 
dynamic Center 


% Increase  in 
Static  Margin 


*(7.92) 

3.59 


*(15.08) 

7.08 


2 3 

(a=128. 59osf)  (o=181 . 69psf ) 


*(15.79) 

11.60 


*(13.39) 

9.09 


♦rigid  model  changes 


Table  28 

Winglet  Effect  on  C^  and  Clq  Derivatives 


Cruise  Condition 
(static  elastic 
models) 




1 

2 

(q=261.67psf) 

(q=128 . 59osf ) 

*(3.25) 

*(3.41) 

.68 

1.76 

*(18.64) 

*(19.56) 

8.00 

13.98 

% Increase  in  Cl 
(more  positive)  a 


% Decrease  in  Cm 
(more  negative)  a 


♦rigid  model  changes 


*(2.89) 

1.30 


*(16.30) 

10.42 


Also  comparing  the  elastic  effect  on  Cma  to  the  rigid  wing- 
let effect,  the  stabilizing  improvement  for  the  elastic 
model  is  significantly  less. 

In  concluding  this  discussion  of  longitudinal  static 
stability,  it  is  interesting  to  observe  the  elastic  effect 
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on  selected  longitudinal  derivatives  versus  dynamic  pres- 
sures shown  in  Fig.  27.  Note  that  Cma  is  the  most  affected 
by  aeroelastic  effects.  Cm^  with  winglets  appears  to  remain 
the  same  as  the  rigid  value  over  the  interval  of  dynamic 
pressures  shown.  No  explanation  for  this  behavior  is  avail- 
able. 

Lateral  Static  Stability.  Table  29  shows  the  winglet 
effect  on  selected  lateral  stability  derivatives.  Note  that 
the  winglets  have  a stabilizing  effect  on  all  these  values. 
In  comparing  the  winglet  induced  changes  with  the  rigid 
changes  that  are  observed,  it  is  seen  that  the  winglet 
effect  on  Cy^  and  Cnr  is  nearly  the  same.  More  changes  are 
observed  due  to  elastic  effects  on  Cn^,  C2,g  and  C£p. 

Fig.  28  shows  the  elastic  to  rigid  ratios  calculated 
with  and  without  winglets  for  the  three  cruise  conditions 
observed.  Ncte  that  for  the  model  without  winglets,  Cj,g 
irs  to  imp.  \..ve  its  ratio  with  the  rigid  value  with 
increasing  dynamic  pressure.  The  explanation  for  this  is 
probably  due  to  the  additional  dihedral  of  the  wing  due  to 
bending.  With  winglets  on  the  elastic  model,  however,  the 
elastic  to  rigid  value  of  C&g  decreases  almost  linearly 
with  increasing  dynamic  pressure. 

Summary.  As  far  as  static  stability  is  concerned,  the 
winglets  contributed  stabilizing  effects  to  CLa»and  Cma. 
Laterally,  the  winglets  were  stabilizing  to  all  the  deriva- 
tives considered. 
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Table  29 

Winglct  Effect  on  Lateral  Stability  Derivatives 


Cruise  Condition 
(static  elastic 
models) 

1 

(q=261 . 67psf ) 

2 

(q=123 . 59psf ) 

3 

(q=181 . 69psf ) 

% Decrease  in  CVr 
(more  \ .gative) 

*(15.52) 

14.56 

.*(15.52) 

14.68 

*(14.29) 

13.4 

% Increase  in  Cng 
(more  positive) 

*(11.76) 

6.0 

*(  5.71) 
3.13 

*(  6.67) 
7.69 

% Decrease  in  CaP 
(more  negative) 

*(30.43) 

4.76 

* (28.57) 
21.05 

* (25.00) 
11.11 

% Decrease  in  C£D 
(more  negative) 

*(11.96) 

3.48 

*(11.88) 

6.67 

*(10.75) 

5.29 

% Decrease  in  Cnr 
(more  negative) 

*(  7.34) 
6.56 

*(  9.71) 
8.19 

*(  6.54) 
5.57 

♦rigid  model  values 


Control  Characteristics  for  Steady  State  Flight 

Longitudinal . The  winglets  on  the  elastic  model  cause 

no  effect  on  the  longitudinal  control  derivatives  Cj^  , 

s 

Cd$s  and  Cm^s  as  Table  III  in  Appendix  E shows.  However, 
since  is  more  negative,  there  is  a requirement  for 
more  negative  stabilizer  trim.  On  the  elastic  model,  the 
winglet  effect  on  is  less.  Consequently,  the  largest 
increase  in  stabilizer  trim  angle  required  is  only  12%. 

This  is  a significant  change  fre  i the  19  to  52%  increase 

S3 
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in  trim  angle  required  by  the  rigid  winglet  model.  Because 
the  trim  requirements  are  less  for  the  elastic  winglet 
model,  there  is  assumed  to  be  little  or  no  noticeable  de- 
crease in  the  range  of  useful  trim  angles  of  attack  and  lift 
coefficients.  The  conclusion  is  that  the  elastic  winglet 
model  experiences  almost  no  trim  penalty.  The  winglet  trim 
effect  is  shown  in  Table  30. 

Using  the  same  equations  discussed  in  part  III,  the 

elastic  model  with  winglets  has  a 9 to  almost  12%  increase 

in  stick  speed  stability,  14^-,  also  shown  in  Table  30. 

dvp 

Again,  this  comes  about,  directly,  from  the  increase  in 
static  margin  caused  by  the  winglets. 

Lateral  Control.  Consistent  with  the  rigid  analysis, 
the  winglet  had  no  effect  on  the  rudder  control  derivatives 
CyfiR,  C£gR  and  Cn5R.  This  is  shown  in  Table  III  of  Appendix 
E. 
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However,  unlike  the  rigid  analysis  where  small  reduc- 
tions in  inboard  aileron  control  authority  were  observed, 
the  elastic  winglet  model  shows  no  change  in  Cy^,  C£gA 
or  Cn§A  comparec  to  the  baseline  elastic  model. 

Steady  State  Maneuvering  Flight.  As  with  the  rigid 
analysis,  all  the  elastic  flight  conditions  studied  here 
are  for  symmetrical  flight  only.  Consequently,  the  maneu- 
vering characteristics  are  confined  to  pull-ups. 

Table  31  shows  the  winglet  effect  in  shifting  the 
maneuver  point  aft.  These  figures  are  approximately  2 to 
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(q=261.67psf)  | (q=128.59osf) 


6.87 


12.03 


9.89 


11.71 


Table  31 
Maneuvering  F: 
State  Pull-up 


1 

(q=261.67psf) 

2 

(q=128.59psf) 

*(7.66) 

3.45 

*(8.26) 

5.88 

*(13.83) 

6.43 

*(15.14) 

11.05 

*(2.89) 

1.48 

* (3.10) 
2.10 

*(15.39) 

*(16.82) 

6.64 

11.89 

4%  less  than  those  observed  for  the  rigid  model.  Using  the 
development  in  reference  6,  maneuver  point  for  a symmetrical 
pull-up  is  defined  as 

M.P.  = xac  - CfnqpScg  (From  Ref  6:5.61) 

4V7 

As  indicated  in  Table  31,  Cmq  becomes  more  negative  with 
winglets  and  from  Table  27  there  is  a winglet  induced  aft 
shift  in  aerodynamic  center  location,  xac.  The  result  from 
the  above  equation  verifies  what  is  shown  in  Table  31. 

The  elastic  winglet  model  has  less  elevator  per  g, 
however,  which  is  also  shown  in  Table  31.  Equations  showing 
how  this  value  is  determined  are  in  reference  6,  page  5.6. 

Summary.  The  control  effects  of  the  winglets  on  the 
elastic  model  are  less  than  the  rigid  model  and  in  some 
cases  negligible. 

Longitudinally,  no  trim  penalty  effect  due  to  winglets 
was  assumed  in  view  of  the  small  changes  in  stabilizer  trim 
angle  required. 

Laterally,  control  effectiveness  in  both  rudder  and 
aileron  was  unaffected  by  the  addition  of  winglets  to  the 
model . 

Consistent  with  what  was  found  in  the  rigid  study,  the 
winglets  resulted  in  increases  in  maneuver  margin  and  in 
elevator  required  per  g for  a symmetrical  steady  state  pull-up. 


The  Winglet  F.f  feet  on  Longitudinal  and  Lateral  Dynamic 
Stability 

Introduction.  Due  to  time  considerations,  this  dyna- 
mic analysis  is  confined  to  studying  just  the  winglet 
effect  on  the  primary  dynamic  aircraft  modes.  Responses 
to  control  inputs  are  not  provided. 

Longitudinal.  The  elastic  winglet  model,  as  shown  in 
Table  32  proves  to  be  less  stable  than  the  baseline  model 
in  the  short  period  longitudinal  mode  for  cruise  conditions 
1 and  2.  Although  cruise  condition  3 is  not  strictly  less 
stable  in  the  short  period  mode,  it  is,  nonetheless,  more 
oscillatory  and  less  damped  than  the  baseline  configured 
model.  In  general,  the  winglet  effect  on  the  short  period 
mode  is  small.  Unlike  the  rigid  analysis  where  the  wing- 
let effect  on  the  phugoid  mode  was  found  to  be  negligible, 
the  elastic  winglet  model  shows  noticeable  improvements  in 
damping  ratios  and  in  times  to  half  amplitude  for  all  three 
cruise  conditions. 

Lateral.  The  winglet  effect  on  all  of  the  lateral 
dynamic  modes  is  stabilizing.  However,  these  improvements 
are  less  than  those  noted  in  the  rigid  analysis. 

The  spiral  and  rolling  modes  are  only  slightly  im- 
proved with  winglets.  But,  like  the  rigid  dynamic  analysis, 
the  dutch  roll  mode  characteristics  are  considered  the  most 
important.  These  are  shown  in  Table  33.  This  winglet 
effect  on  dutch  roll  stability  comes  primarily  from  winglet 


Table  32 

Winglet  Effect  on  Longitudinal  Dynamic  Modes 


Cruise  Condition  j 

1 1 

2 

3 

(static  elastic  models) 

M = .77 

M = .77 

M = .422 

i 30,000  ft  1 

1 45,000  ft 

10,000  ft 

I 

Short  Period 


% Increase  in  Un- 
damped Natural 
Frequency , wn 


3.47 


4.45 


% Decrease  in  Damp- 
ing Ratio,  cs.„. 


3.89 


6.25 


4.11 


% Increase  in 
Time  to  Half 
Amplitude,  T 1/2 

% Decrease  in 
Period,  T 


.49 

.45 

0.0 

4.34 

6.57 

5.79 

Phugoid 

% Increase  in 
Undamped  Natural 
Frequency , wn 

% Increase  in  Damp- 
ing  Ratio,  gp 

% Decrease  in 
Time  to  Half 
Amplitude,  T 1/2 

% Decrease  in 
Period.  T 


4.33 

1.45 

3.00 

6.60 

1.69 

4.08 

10.07 

3.15 

6.35 

4.10 

1.57 

2.68 

Cruise  Condition 
(static  elastic  i 
Dutch  Roll 


% Increase  m Un 
damped  Natural 
Frequency , “n 


% Increase  in  Damp 
ing  Ratio,  5d.R. 


% Decrease  m 
Time  to  Half 
Amplitude,  T 1/2 


% Decrease  in 
Period,  T 


1 

M = .77 
30,000  ft 

2 

M = .77 
45,000  ft 

3.75 

2.02 

11.46 

5.48 

13.05 

7.51 

3.5 

1.90 

W 
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V.  Conclusion  and  Recommendations 

Using  geometry  and  aerodynamic  data  on  the  KC-135  and 
the  Boeing  winglet,  a mathematical  rigid  and  elastic  FLEX- 
STAB  model  was  successfully  constructed.  The  correlation 
between  the  baseline  model  and  available  flight  test  data 
was  considered  close  enough  to  regard  the  model  as  a good 
approximation  of  the  actual  aircraft.  This  flight  test 
comparison  was  made  for  both  the  elastic  and  rigid  models 
where  the  elastic  structural  definition  consisted  of  a 
modified  707-320B  model.  The  winglet  geometry  was  modeled 
in  detail  with  both  thickness  and  twist  definitions , 
however,  maintained  as  a rigid  structural  body. 

In  the  rigid  study  of  the  winglet  effect,  the  baseline 
model  was  compared  to  the  winglet  configured  version  for 
five  flight  conditions  which  represent  the  KC-135  mission 
profile.  The  winglet  effect  on  performance  for  the  rigid 
model  showed  a 2.32  to  7.89%  decrease  in  total  airplane 
drag  with  an  accompanying  2.47  to  8.58%  increase  in  lift 
over  drag  ratios.  However,  an  increase  in  negative  stabil- 
izer trim  angle  was  consistently  noted  with  winglets.  Also, 
the  range  of  useful  trim  lift  coefficients  and  trim  angles 
of  attack  was  reduced.  This  result  was  directly  due  to 
increases  in  CLa  and  to  the  winglet  induced  aft  shift  in 
the  aerodynamic  center  of  the  wing,  which  resulted  in  a 
15.85  to  28%  increase  in  static  margin.  Laterally,  the 
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winglets  were  stabilizing  to  Cn^.  In  aircraft  control,  the 


winglets  had  no  effect  on  elevator  authority,  but  some 


small  reductions  in  inboard  aileron  effectiveness  were 


noted.  The  rudder  was  unaffected  by  adding  winglets  to 


the  model. 


Dynamically,  the  rigid  winglet  model  displayed  a 6 to 


almost  12%  reduction  in  the  damping  ratio  of  the  short 


period  mode  along  with  8 to  13%  increases  xn  undamped  natu- 


ral frequencies.  No  effect  due  to  winglets  was  seen  on  the 


phugoid.  Laterally,  the  rolling,  spiral  and  dutch  roll  modes 


were  all  more  stable  with  winglets.  In  particular,  de- 


creases in  spiral  divergence  were  noted  along  with  a 


10%  reduction  in  the  rolling  mode  time  to  half  amplitude. 


The  dutch  roll  winglet  effect  was  regarded  as  the  most  sig- 


nificant, however,  with  3 to  12%  increases  in  damping  ratios 


and  5 to  14%  decreases  in  times  to  half  amplitude. 


Also  included  in  the  rigid  analysis  was  an  investiga- 


tion of  the  lateral  stability  characteristics  of  the  winglet 


model  with  a truncated  vertical  stabilizer.  This  resulted 


in  less  stable  side  force  and  yawing  moment  derivatives 


than  for  the  baseline  model.  Reductions  in  rudder  authority 


below  the  baseline  configuration  were  also  noted.  Dynami- 


cally, the  short  tail  model  was  less  stable  in  dutch  roll 


over  the  full  tail  model  without  winglets.  No  effect  on 


the  rolling  mode  was  seen.  The  most  remarkable  result  was 


the  large  increase  in  spiral  stability  seen  with  the  winglet 


short  tail  model  over  both  of  the-  full  tail  conf iqurations 
compared. 

Elastically,  winglet  performance  improvements  were 
degraded  primarily  by  the  increase  in  "washout"  at  the  wing 
tip.  Here,  winglet  drag  reductions  for  the  three  cruise 
flight  conditions  studied  amounted  to  only  1.3%,  with  only 


1 

I 


I 

£ 

I 


a 1.3  to  1.54%  increase  in  L/D.  Trim  requirements,  however, 
were  much  less  than  those  found  in  the  rigid  analysis.  Con- 
sequently, the  elastic  winglet  model  showed  no  noticeable 
trim  penalty.  The  winglets  were  stabilizing  longitudinally 
with  7 to  11.6%  increases  noted  in  static  margin.  Laterally, 
the  elastic  winglet  model  revealed  the  same  general  trend 
in  improving  stability. 

Dynamically,  the  winglets  had  less  of  an  effect  in 
reducing  the  short  period  damping  ratio  on  the  elastic  model. 
Winglets  also  had  a stabilizing  effect  on  the  phugoid  mode 
with  small  increases  in  damping  ratios  and  3 to  10%  decreases 
in  times  to  half  amplitude.  Lateral  dynamic  improvements 
due  to  winglets  were  not  seen  in  the  spiral  or  rolling 
modes.  However,  as  with  the  rigid  analysis,  the  dynamic 
characteristics  of  the  dutch  roll  mode  improved  with  wing- 
lets. Here,  5.48  to  11.46%  increases  in  damping  ratios 
along  with  7 to  13%  decreases  in  times  to  half  amplitude 
were  seen. 


ij 


a 


Generally,  the  overall  winglet  effect  on  performance 
for  the  models  used  in  this  study  was  not  as  great  as  that 
achieved  in  wind  tunnel  tests.  Nevertheless,  stability 
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Airplane  Small  Perturb?. ‘ ion  Equations 
Used  in  Dvnamic  Time  Ilistorv  Analysis 


The  dynamics  of  this  report  are  concerned  only  with 
small  perturbations  relative  to  a steady  state  flight 
condition  in  which 

a.  no  initial  bank  angle  exists  ( ^ = 0)  , 

b.  no  initial  side  velocity  exists  (V^  = 0) , and 

c.  no  initial  angular  velocities  exist: 


£*1  = Ql  = Rl  = 4*1  = 6i  = 4>i  = 0 


(From  Ref  6:2.36) 


Therefore  the  equations  of  motion  become  longitudinally , 


m(u  + W^q)  = - mgScosO^  4-  fa  + fT 

X X 

m(w  - ulQ)  = “ mg0sin8^  + fAz  + fT 


q = 0 


(From  Ref  6:2.36) 


and  laterally. 


m(v  + Ujr  - W^p)  = 

zxxP  ’ xx9P  = 
- Z„  p = 


mg^cos©!  + fAy  + fTy 


+ 


nA  + nT 


= $ - r tan8^  (From  Ref  6:2.36) 


These  equations  are  then  written  in  the  stability  axes 
system.  Therefore,  Wj_  = 0 in  both  equation  sets  (1) 
and  (2)  . 

The  inertial  transformation  from  an  arbitrary  body 
fixed  reference  axis  system  to  the  stability  axes  system 
is  accomplished  using  this  matrix  equation: 


Note  that  I.,.,,  is  invariant  when  a coordinate  trans forma- 
tion  consists  of  just  a rotation  about  the  airplane 
y-axis  (Ref  6.5.43). 

From  the  derivation  of  the  perturbed  state  forces  and 
moments  given  in  section  4.2  of  chapter  4 in  part  one  of 
Reference  6,  the  complete  small  perturbation  equations  of 
motion  for  steady  state  flight  become  those  given  in  Table 


I.  Note  that  these  equations  assume  that  CD*,  Cn  and 

a ' uq 

(Cm  + 2 Cm  ) are  approximately  equal  to  zero  for  the 
2u  Z1 

longitudinal  set  and  that  cyy  cZg  and  cnjj  are  approximately 
zero  for  the  lateral  set  (Ref  6:4.113).  In  the  FLEXSTAB 


characteristic  equation  rooting  analysis,  however,  these 
derivatives  are  not  assumed  to  be  zero.  This  accounts  for 
the  small  differences  noted  between  the  eigenvalues  of 
FLEXSTAB  and  those  calculated  using  CALMAT. 
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In  Roskam's  discussion  of  dynamic  stability  and 
response,  he  develops  what  he  calls  "dimensional  stabil- 
ity derivatives."  These  allow  simpler  forms  of  the 
equations  in  Table  I.  The  longitudinal  equations  are 
divided  by  m or  IVy  where  applicable  and  the  lateral 
equations  are  divided  by  m,  Ixx  and  1Z2»  respectively. 

This  results  in  the  following  forms  of  the  perturbation 
equations  in  terms  of  "dimensional  stability  derivatives": 
Longitudinal — 


u = - g0cosG1  + XUU  + XpuU  + Xaa  + X^e6e 
w - Uiq  = - gesinGiL  + Zuu  + Zaa  + Z*a  + Zqq  + Z6e6e  (4) 

q = r^u  + M^uu  + Maa  + ^Taa  + + M^q  + 

§ = q 

(From  Ref  6:6.19) 

Lateral — 


v + Uxr  = g4>cos0^  + Ygg  + Ypp  + Yrr  + Y6a5a  + Y6r6r 
p - A^r  = Lg8  + Lpp  + Lrr  + L5a6a  + l<5r6r 
r - Bxp  = n6B  + NT(JS  + NpP  + Nrr  f N6a6a  + N6r5r 
$ = p + rtan6]_ 

where  = IX2  - Ixx,  = IX2  t i22 

(From  Ref  6:6.47) 


The  dimensional  stability  derivatives  of  equations  (4)  and 
(5)  are  defined  in  Tables  II  and  III  of  this  appendix. 

Equation  sets  (4)  and  (5)  can  now  be  arranged  in 
first  order  matrix  form  as  shown  in  Table  IV.  In  this 
state  space  form 


107 


W 

I 

£ 


Table  II 

Longitudinal  Dimensional  Stability 
Derivatives 


9 

if 

5 

5 

! 

i 

I 


- .**W*,f  VW>*a»,  ;,*&»**«■« 


Table  IIJ 

Lateral-Directional  Dimensional 
Stability  Derivatives 


_qiSCyB  (ft  sec-2) 


qiSbCyp  (ft  sec-1) 
YP  = 2mUl  “ 


_ «Sb^.  (ft  sec-1,  _ (ft s-^2d--l, 

2mUi  6A  m 


qiSCv5  (ft  sec-2  or_ 
y6R  = — L— i-  ft  sec"2  deg 


(ft  sec-2  or  i ^ q^SbCjtp  (sec  ) 

ft  sec-2  deg-1)  Lg  = 


qiSb2Cap  (sec-1) 

= “2i^xUr 


qiSb2C*r  (sec  -1) 
Lr  ~ 2lxxuY” 


Ne  = -r^" 

qiSb2Cnp  (sec"1) 

^O  " ot n. 


(sec  2 or 
sec-2  deg- 

LSr  = 

qiSbCjigR 

^•xx 

(sec-2) 

air 

c^SbC^ 

Izz 

vsec 

sec-2  deg-1) 
(sec-2) 


_ qiSb2Cnp  (sec"1)  qi Sb2Cn_  (sec-1) 

P ' 2IzzUi  Nr  = 2IzzUi 


qjSbCn^  (sec"2  or  sec"2  deg"1) 

NSA  = 


q^ShC^  (sec-2  or  sec-2  deg"1! 


NfiR  I 


(From  Ref  6:6.46) 


**  jus*  «- 


Table  IV 

Aircraft  Equations  of  Motion  in  First  Order  Matrix  Form 


Longitudinal — 


u 

XU' 

xa 

0 

-geos  6]^ 

• 

Zu 

Za 

Ul+Zq 

gsin01 

a 

U1-Z& 

Ul-Za 

Ul-ZJ 

Ul-Za 

• 

q 

* , MaZu 
Mu  + Ul-Z^ 

MaZa 
Ma  + Ul-Zi 

„ M;<Za+  Di) 

q Ul-Za 

M&gsin0 

ui-z& 

0 

0 

1 

0 

— X6E 

X$t 

Z5e 

Zfiip 

Ul-Za 

Ul-Z& 

MSv  + Maz6E 
Ul-Za 

MaZST 

M5t  + Ul-Za 

__  0 

0 _ 

Lateral — 


• 

A 

Y6 

Yr  _ i 

geos  01 

6 

P 

Ul 

Ul 

Ul 

Ul 

• 

ainb+lb 

AlNp+Lp 

AiNr+Lr 

0 

p 

p 

1-Alfll 

1-A1B1 

1-AlBi 

• 

3 

bilb+nb 

BlLo+Np 

BlLr+Nr 

0 

r 

JL 

1-AlBl 

1-AlBl 

1-AlBl 

• 

* 

0 

1 

tan  6i 

0 _ 

- <J>_ 

Y5a  y$R 

l”5a  l^R 

5A 

N|a  NjR 

5r 

0 0 


where 

Al  = tX2  and  Bl  = lb'- 
Ixx  Izz 
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Appendix  B 
Program  CALMAT 

CALMAT  is  a digital  computer  program  that  computes 
and  prints  the  stability  and  control  derivatives  matrices 
for  the  longitudinal  and  lateral  matrix  equations  of 
Table  IV  in  Appendix  A.  The  program  also  uses  A.F.  Flight 
Dynamics  Laboratory  EISPACK  subroutines  to  calculate  the 
longitudinal  and  lateral  eigenvalues  and  respective 
eigenvectors  of  the  stability  derivatives  matrices.  This 
function  of  EISPACK,  used  by  CALMAT,  is  in  effect  a cha- 
racteristic equation  rooting  analysis.  After  the  roots 
are  calculated  by  EISPACK,  CALMAT  then  computes  and  prints 
the  undamped  natural  frequency,  damping  ratio,  period  and 
time  to  one  half  amplitude  of  each  mode.  Any  unstable 
roots  are  also  noted  and  the  time  to  twice  amplitude  is 
computed  and  printed  for  these. 

The  program  can  be  used  to  calculate  either  body 
axes  data  matrices  or  stability  axes  matrices.  As  used 
in  this  report,  all  data  input  to  CALMAT  are  in  the 
stability  axes  system. 

Data  required  for  program  execution  makes  up  eleven 
cards  specified  in  Table  I of  this  Appendix. 

The  matrices,  eigenvalues  and  eigenvectors  computed 
by  CALMAT  are  used  in  computing  the  modal  time  history 
and  time  history  due  to  controller  inputs  of  each  of 
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the  longitudinal  and  lateral  state  variables. 

The  program  listing  is  shown  in  Table  II  and  the 
program  flow  chart  in  Table  III. 


Table  I 

CALMAT  Data  Card  Format 


Card  nu. 

Data 

1 

9/  OB/  Ti 

2 

q,  Ui,  Wl 

3 

m,  S,  c,  b,  dT 

*4 

!xx»  !yy»  Jzz 

5 

cLl#  cLct'  CLa'  CLq* 

CL«e.  Cl U 9 

6 

cDl»  CD(X,  CDu,  CDfie 

7 

^-mi  > Cma  , , Cniq , 

crau/  cm5e 

8 

cmTa 

9 

^Ya'  cYp'  CYr f CySA' 

CY6R 

10 

cIq>  Cip,  c5.Sa' 

cUr 

11 

cnp/  CnTo / Cnp,  Cnr, 

cnfiA'  Cn6R 

Dimensions 

radians,  pounds 
psf,  fps 

slugs,  ft2,  feet 

slugs-ft2 

rad--*- , sec-rad-1 

rad-1,  sec-rad-1 
rad*"1  f sec-rad-1 


rad-1,  sec-rad-1 


rad-1,  sec-rad-1 
rad-1,  sec-rad-1 


♦Mass  moments  of  inertia  must  be  input  in  the  principal 
axes  system  of  the  aircraft.  The  program  uses  the  trans- 
formation of  equation  (3)  in  Appendix  A to  compute 
inertia  terms  in  appropriate  axes  system. 
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Table  II 

CALF-1AT  Program  Listing 

PR0'.7AH  CAL«<VT(P!POT  = /ia  .D'JT-’JT  , 24  3S5  = I N3  JT  , f 43£&  = 3'JT0,U) 

oih'msion  p(<*, ?i , cas r < «•! 

re  a-.  m,  t zT,  ixz,“-j,  -» t 'j,  - a ,Mri  ,Mn,-a,'Oi  ,nT,.  i,l3,'.  j , 

lL9A,0?.,‘*T»NrT,'43,‘l3,-|-:;,'nr,A  (*,«)<  om£gm  (4> , ze  tau)  , th4lf<  ♦> , 
1^:131, °.CIO?,  =£3tOn(  M , P3:%  T’C.1 

COMPLEX  Z1WV.  (*)  ,i:o  i:r.C.,  ~),.ORK<4>  .EIG.CEIGVAL 
reai  a,  (oaseo 

I-  (*nF{5Li'<°jn  .*:•  .i)  sn  = ••''D'lPurATiC'j;  ;oh°l£T£" 

FOP.^AT  1 6 A 1 0 ) 

i:ojht=o 

READ*,  TH£TA,Al.3HA,Tl 
REAn*,D,t|,H 
READ»,H,3  ,C«AD,n,9T 
REA3*,  ixx,  :yy,  127 
AHG'-ES,  03  *•$<■  A_3HA 

SIXX=GDS(  AUGUE)  ‘ *2*  IXX*9r‘}<4'43L£)*»2*T!7 
SIZZ=SIN(  ANGLE)  **2*IXX*CDS<  A'iGLE)  **2*122 
IX7=.5*SIM(?.M4NGl£) )*(IXX-I72) 

IXX=3IXX 

IZZ=SIZZ 

****** LONGITUOI'JAL  CALCULATIONS 
Cl=D*S/H*U) 

C2=(o.-i.)*:i 

C3*02*U 

C%=02*CMAC/2. 

C5=Cl*M*CMAC/lYr 

■ CS=0>*U 
cz*:5*chac/2. 

■ .R£AD*,CL1  »CL  4»CLAD,  CLDfC.DEtG.’J 
REA9*,C9i,COA,CVJ,CDO£ 

^;A3*,CM1  ,C1  A.GHAO,  C'n,:-4U,3'fD£ 

R£A3*,CHTA 

Srxiscot- 1**32.174)  •sr'HTHirAj/cn*?) 
crx'j=(o.-2,)  *:txi 
01T1=  <Q.-  1«)  *ZU 
mtu=(0.-i.)  »:rxu*or/ci4C 
CU=02*(CTJ+2.»CD1> 

XTU=Cl*  (CTXU  +2,  »CTX  J.) 

-xa=:3»{coa-;ld 
X3E*C3*C09E 
Z’J=C2*<CLU*2.*CL1J 
za*C3*  (cla*:oi) 

ZA0=C4*CL  40 
ZQ*C4*CU1 
Z3EO3*CL0E 
>u*:s*cc'ui+2.*C'U> 

MTU=05* (3  MTU *2.*CHT1) 

HA*:3*CM4 
iZA=:6*CMTA 
XA0=37*C*1  AO 
Ma=07*CH1 
«3E=:6*C10E 
OT1=T1/10  00. 

COOT* (0,-1.) *:oi/OTl 
XOT=C3*COOT 
CMOTsCMTl/OTl 
M0T*C6*CM0T 


I 


TaMe  II  (cont.) 


a(i,ii=xu»xtu 

A(1,?1=XA 

««i.3»=C0.-l.)*W 

A«1,,4)=Ci).-1.)  M?.l?4»033<  TH-TM 
OiH=J-r AO 
At?, 1 >=ZU/0E M. 

At  2,21  = 71/0" M 
At2,31=  ('J*?01/0£N 

At2»4)  = (3.-1,)M2.174*3IN(J''<“7A1/3EN 

A{3,1)=MU4-mTU»-1A0*ZL'/0£N 

A (3,  2)  = MTA+ha*-  1A 0 *7  A/  0£N 

At  3,  3)  = HT4MA  O*(ZO*iJ}/0£N 

A(3»4)=tO.-l.l*'lA‘)*32.t74*SI‘l(TH£TAl/DEN 

A 1 4, 1 ) = 0,  0 

A 14,21  = 0. 0 

A 1 4,  3)  = 1.  0 

At  4,  41  = 0,  0 

8(l,ll=XOE 

9tl,2)=X0T 

Bt2,ll=20E/0EN 

Bt  2,21  = 0.  0 

B<  3,  1 1 = M3£+m  0»7Q£/  0£N 

B t 3, 2 1 = H0  T 

Bt  4, 1 1 = 0. 0 

8(4,21=0. 0 

PRINT* 

‘PRINT* 

PRINT* 

°RINT*,“ 

PRINT  7,  tCASEdl  , 1=1,61 

7 FORMAT  tllX,  6A101 

Pm*," 

PRINT* 

PRINT* 

PRINT*,  “BASS  MOMENTS  OF  INERTIA  I >4  STABILITY  AXIS  3rST£M“ 
PRINT* 

• PRINT*,  “IXX=“,IXX,“  IZZ=“,IZZ,“  IXZ=“,IXZ 
PRINT* 

PRINT* 

PR INT*, "LONG ITU3 INAL  A MATRIX" 

PRINT* 

00  20  1=1,4 

print  s,  t At i, ji  , j=i,4> 

8 FORMAT  (5X,4F10.41 
PRINT* 


ao . continue 

PRINT* 

PRINT*,  "LONGIT'JOINAL  9 MATRIX- 
PRINT* 

03  30  1=1,4 

PRINT  9,  tBt I, Jl  , J = l,2) 

. % FORMAT  t5X,  2F10.41 

PRINT* 

30  CONTINUE 

CO  TO  103  U6 
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I 


1 


I 


| 

I 

I I 


i 

i 
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Table  II  (cant.) 


c.,.*,,»..;_arr^flL  CALCULATIONS 
2 C8=")»S/m 

C9=C6*?9/  (Z.  *'J) 
0100*3*93/1 XX 

cno  in  ^/t  ?..  * j) 

Cl ?0*S*  3 9/1 71 
C13=D12*39/(  2.*\1) 
E=I*Z/Ixx 
p*IXZ/I7Z 


READ*, 

CY 

0, 

CY°,  CY  R, 

CY  DA , 

CYDR 

RE  AD*  , 

Cl 

3, 

C1P,  C1R, 

C1DA, 

01DR 

READ*, 

CM 

9, 

C'lT  3 * CNP, 

CNR, 

CMDA , CNOR 

Y3  = D8*CO 
yp=:9»cvp 
YR=C9*CYR 
YDAOO»CY  DA 
YDR=C8*  CY  OP 
L3=C10*C13 
L9=C11*C1 ° 

LR=C11*C1R 
L9A=010*C10A 
LDR=C10*C10R 
N3=S12*CN9 
NT  9= D 12 ‘SNTP 
NP=C13»CNP 
NR=C13«CNR 
'N>A=C12*Ct<0A 
NDR=C12*CN0R 
H=i.-E*F 
Mil  11 =Y9/U  . 

A(1,21=IYP+W)/'J 
A(l,31=YR/'J-t.n 
»a,4)=32.17-»CDS<THrTA)/U 
A(2,11=<E* <N-3*Nr9)*L9)/H 
A(2,21=<E*NP+lP)  /H 
A(2,3)=(E*NP+LR)/m 
A(2,41  = 0.0 

A(3,11=N3>NT3*F*  <£*  (N3*NT3>  *L9)  M 
A(3,21=NP*F*  <E*N°+LP)  / H 
A(3,31=NR*F* <E*NR+LP> /H 
A(3,41=0, 0 
A(4, 1>  =0. 3 
A<4,21  = 1.0 

A (4, 31=SIN(THITA1/C0S (THETA) 
A(4,41  = 0.0 
'Btl.l)=Y0A/U 
3(1, 21 =Y0R/U 
3(2,ll  = (E*NnA«-L9A)/H 
8(2»21=(E*N0R*LDR)/H 
3(3, 1)=NDA+F*9(2, 1) 
8(3,?1=N0R+F*3<2,2) 

3(4,11  = 0.  0 
3(4,21=0,0 

HINT* 

PRINT*,  -LATERAL  A MATRIX- 
'PftlNT*  • U7 


Table  II  (cont.) 


10 

40 


it 

45 


00  I =t  1 4 

p^x-if  to*  (4 
F0P.Jtt”  C 5 < » 4' n • 

P3  t-J  T * 

-ourt'MOE 

ss;:,  -uat^v.  s -mx*- 

00  +5  I =1 > 4 , , ,, 

P^IMT  Ilf  (q  (I» J)  ’ J'i’2 
P0P.-14T  (5X»  2r  to « -> 

97.V  IT* 

e5W”«Swl^»-US  GO^umiOSS 
103  X:OJMTstCOUMT*l 
Hs  4 
IHOsl 

• SST/.,  -cwwcrsnsTis  swmoH  ootug  *«t»SXS- 

2°ITE(5^301) 

FOR^flTdXj&t  1P016.S)) 

"*  g£'Sr igW,».».™.™«-''£II''K,"0RO 

W?lTs(6>3)  (*IG-/4t(I)  fl-1*  4) 
p«*»T  ieo,<i»t*t»^  -,23X),'> 

180  F0R4AT{//,T13,'..<  -IGV"  *“* 

03  5 1=1. K t-1  M) 

Nirs(6«&>  (n3'/£cii,J>»J-i»S) 

6 F3<mT(lX,9<l:>Si6*d)) 

* CONTINUE 

03  400  K=1»N 
EIG=EIGVAL(X) 

^F^tRSIGl.GT.O.O)  GO  T0  399 

SSSSSSilS..* 

g*?isssts”!“j«»«"  ” !,j 

EtG=nGVAUK)*»iIGVAL 

PVPsSI6 

OHEG'<(K)  = SORT(5pOr)> 

ZiTA<K)  t0«GM  «)  -Son  Cl,-XSTa«l  **2>  > 

SSfiKl*  ,693dOMEGMU)*?*-TA«l) 

P*INT* 


loO 
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170 


IfXETUK) 

nSTrN-T<-,<.->  »-SE=S-,- 

U«,-SECS- 
P*XNT* 


Z£TA'  *<»~) 


11. 


^ faiiS'Kja  a irastc.-* 


Table  III 
CALMAT  Flow  Chart 


© 


CALCULATE 
Along  & Blong 


m,S,c,b,dT 


Ixx'Iyy'Izz 


CALCULATE 

Ixzs# 

I2ZS 


PRINT 
Headings  & 
Along  & Blonc 


CALCULATE 
Eigenv  alues 
& 


CALCULATE 
Aiat & Blat 


PRINT 

PRINT 

Eigenvalues  & 

Headings  & 

Eigenvectors 

L 

Alat  & Blat 

g> 
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Appendix  C 


Program  STATVAR 

STATVAR  is  a digital  program  that  numerically  inte- 
grates the  longitudinal  and  lateral  equations  of  motion 
(first  order  form.  Appendix  A)  of  an  arbitrary  aircraft, 
producing  the  values  of  the  state  variables  as  a function 
of  time.  Inputs  required  are  the  integration  time  step, 
the  maximum  time  value  desired,  the  dimensional  stability 
derivatives  matrix,  the  control  derivatives  matrix  and  the 
initial  values  of  the  state  variables.  Different  schedules 
and  values  of  control  inputs  may  be  inserted  by  means  of 
two  control  function  subprograms.  Output  consists  of  a 
time  listing  of  the  state  variables  as  well  as  CALCOMP 
plots  of  the  longitudinal  and  lateral  time  response.  The 
value  of  the  time  step  T may  be  varied  within  a limited 
range  without  seriously  deteriorating  the  accuracy  of  cal- 
culation, by  adjusting  the  number  of  terms  taken  in  the 
expansion  of  eAt. 

Using  the  first  order  form  of  the  equations  of 
motion  as  developed  in  Appendix  A,  they  can  be  written  as 

X = A xT  + B u (1) 

where  A and  B are  the  matrices  of  stability  derivatives 
and  control  derivatives  respectively,  x is  the  column 
vector  of  the  state  variables  and  u is  the  column  vector 
of  the  control  variables. 
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Taking  the  Laplace  transform  of  equation  (1)  with 
zero  initial  conditions  and  rearranging: 


sX(s)  - AX (s ) = BU (s) 


X(s)  = (si  - A)  _1BU  (s) 


(2) 


Finding  the  homogeneous  solution  of  equation  (2) 


X(s)  = (si  - A)“1Xf 


and  taking  the  inverse  Laplace  transform  yields  equation  (3) 


:(t)  = L-1[(sI  - A)-1]x0  = eAt 


'x0 


(3) 


where  eAt  = L-1(sI  - A)~^.  Finding  the  particular  solution 
of  equation  (2)  for  an  impulse 


X (s) 


imp 


= (si  - A) -1B 


{:} 


Let  X(s)imp  = H (s) * , where  H(s)*  is  a matrix  of  transfer 
functions  between  each  of  the  state  variables  and  control 
variables : 


e.g. 


H (s) * = (si  - A) -1B  = 


G 
.e 


u<5<»  uU6t 


G®6e  G06t 


Taking  the  inverse  Laplace  transform  of  H(s)*  yields 
equation  ( 4 ) : 

x(t)imp  = L_1[(sT  - A)"1]!  = eXtB 
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(4) 


a 


*-  - <•'<*«-'**?.  v*f1»ywr*vj- 


It  can  be  shown  that  an  arbitrary  control  function  can 
be  represented  by  a summation  of  impulse  functions  over 
the  time  interval  in  question.  Therefore,  the  particular 
solution  for  an  arbitrary  control  input  can  be  written  as 


t 


x(t)p  = £eA(t  “ t^Bu  (T)dT 


Summing  the  homogeneous  and  the  particular  solution  gives 
the  complete  solution: 


t 7, 


x(t)  = eAtXQ  + oeA^  ~ t*Bu  (x)dT 


(5) 


It  is  this  equation  that  is  numerically  integrated  by 
STATVAR. 

Selecting  a time  step  T,  the  homogeneous  solutions  at 
times  tg  = 0,  t]_  = T and  t2  * t^  + T are 


x(tg)  = X0, 


x(ti)  = e xg,  and 


:(t2)  = + T^xg  = eATx(tx 


This  implies  that  x(tn)  = eATx(tn-i) 


or 


x[(k  + 1 J'T]  = eATx(kT) 


For  the  particular  solution 


T - — t — 

x(T)d  = ^eA(T  - t)b  u (x)dx  = eAT*e'ATBu  (T)dr  (6) 
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For  small  values  of  T,  u(x)  can  be  approximated  by  a 


constant  u(0).  Therefore 


x(T)p  = (e^e'MdT}Bu(0) 


The  complete  solution  is  now 


x(T)  = eATx(0)  + (eM^e‘^TdT)Bu  (0)  and 


x{2T)  = e2A>rx{0)  + e2AT  Oe~ATdxB  u 

t 0 


2 T “ __  _ 

+ ~ATdxB  u 
T 


U (T)j 


This  can  be  written  more  simply  as 


x(2T)  = eATx(T)  +^eAT^e  AxdxJ  B u (T) 


or,  in  general 


x[(k  + 1)iT]  = 


eATx(kT)  + (eAT'e"Axdx)B  u (kT) 
0 


x[(k+  1)T]  = eSTxxkT)  + A"1  (eAT  - I)B  u (kT)  (7) 


where  eAT  can  be  calculated  by  a Taylor  series  expansion: 


AT  — — T^  — -j 

e = I + AT  + A 2 =-r  + A3  i-r  + 


Program  STATVAR  carries  out  the  integration  of  the  time 
domain  equation  (5)  by  iterating  solutions  of  equation  (7) . 


•a&sa&s® 


yaL.’ja.w 


The  accuracy  of  the  integrated  solution  of  the  state 
variable  equations  is  affected  by  the  integration  tine 
step  T,  the  number  of  terms  calculated  in  determining  eAT 
from  the  Taylor  series  expansion  and  by  the  type  of 
control  function  used. 

The  error  in  the  calculation  of  eAT  is  determined  by 
the  size  of  the  time  step  and  the  number  of  terms  which 
are  taken.  In  general,  to  maintain  a desired  maximum 
error  as  the  size  of  the  time  step  is  increased,  more  terms 
must  be  taken  in  the  expansion.  For  example,  for  a T = .1 
seconds,  a four  term  expansion  will  yield  an  error  on  the 

order  of  A4  Zl  = A4  x 10“6.  If  the  time  step  is  increased 

4! 

to  T = .5  seconds,  then  a four  term  expansion  will  yield 
an  error  on  the  order  of  A4  x 10“3.  Therefore,  more  terms 
must  be  taken  in  order  to  maintain  the  desired  error  margin. 

Also,  in  developing  the  iteration  of  equation  (7) , it 
was  assumed  that  u(t)  was  essentially  constant  over  the 
time  step  T.  Consequently,  a rapidly  varying  control 
function,  such  c.s  a sinusoidal  input,  will  require  a smal- 
ler time  step  than  would  an  essentially  constant  control 
function,  such  as  a step  input  or  zero  input. 

See  Tables  I,  II  and  III  for  input  requirements,  flow 
chart  and  program  listing. 
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Table  I. 

STATVAR  INPUT  REQUIREMENTS 


I FLAG: 


0: 

T: 

TMAX: 

I UNITS : 

M: 

N: 

MB: 

NB: 

AMAT ( I , J ) 
BMAT ( I , J ) 


0 < indicates  longitudinal  case  only 
0 > indicates  lateral  case  only 

= 0 indicates  both  longitudinal  and  lateral  case 
indicates  number  of  terms  to  be  calculated  for  e 
time  step  (seconds) 
maximum  time  (seconds) 

= 0 indicates  input  angles  in  degrees 
= 1 indicates  input  angles  in  radians 
number  of  rows  of  A matrix 
number  of  columns  of  A matrix 
number  of  rows  of  B matrix 
number  of  columns  of  B matrix 
elements  of  A matrix 
elements  of  B matrix 


AT 


X(I)  : 


*Note:  The  units  of  the  A & B matrices  must  be  compati- 

ble with  the  units  of  X(I),  FU1  and  FU2,  e.g. , 
radians  with  radians  & degrees  with  degrees. 

initial  conditions  of  state  variables 
Control  Function  Subroutines 


Various  control  functions  may  be  inserted  into  the  program  by 
inserting  different  control  function  subroutines. 


FU1: 

FU2: 


elevator  and  aileron  control  functions 
throttle  and  rudder  control  functions 


Some  example  control  subroutines  are  included  with  the  pro- 
gram listing. 
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Table  II 

STATVAR  Flow  Chart 


Table  III 

STATVAS  Program  Listing 
PRO! RAH  STAT'/a?  CIN3UT  = /3C,  1 J"3 JTr/l^Q  ,P(_0T ) 

OIVUGIOM  AT  ( + ,•*)  ,0  lATU.Rl.XCu)  ,UK2)  , AI  C4,V  ),  1.4(4)  , "1  C..)  , 

1 FATC4,  4),  APATC4.4)  ,0X  (4,1)  ,AHU  ,M,H<  4,2)  ,H'J(U,  1>  , X=(u) 

i UP  C2) 

COM“OM/UFU>T  /'i!cUf 

COMMOM/L  NO/ Ti  Cl  05)  ,11(315)  ,«C»n  5),  M f!33>  ,TH(  3Q  *5>  , fl£  C 305)  , OT  ( J05) 

1 

C0N"0N/L4  T /T2CT01;)  ,0(305)  ,3  dOS)  ,9£CT05>  ,3H(T0  5>  ,9W{T05)  ,?R<T05) 
1 L4T  C 2 ) 

INTEGER  0 
N»i  = 0 
NP?  = 0 

t PlftO^IPL'^O.T.T^av.TimTS 
IFroKStmur  > iso,  io 

10  IFCIFLAG.  Lt.  0)  >£40  5 , CLONC  7 > , 1 = 1 , 2 ) 

IF  CTrLAG.  G”.  0)  READ  5,  (LOff  n,I=l,») 

5 FOR'4  AT ( ?4  1 P) 

REA0* ,M,M ,mo  ,>n 
IF (“Or(FL I^’FUT)  > t^n,->o 
20  RE  A'-*,  < (A  hat  (I , J ) , J = 1 , N>  , I = 1 ,") 

REA3*,  ( C3“A~  (I,  J)  ,j=l,*n>  ,T=1,H9> 

NPL9T  = TMfiy/T/Vn. 

NpTS  = TH AX/T/N3LOT 
Ic(NpTS.GT.3<)2>  N3L0T  = NpL3m 
CALL  rATH  AT  C AMAT  ,EAT,1£4T  ,H,N,  0,T) 

03  30  1=1, M 
00  *0  J=l,U 
30  Aid,  JJsA^ATU,  J) 

CALL  MTIN(AI,M,3,U,MA) 
lc«n,--0.0.)  GO  TO  IPO 
REA3*, <XCI),I=1,M) 

TTOT=0. 

IPLOT  = 0 

UTC1)  = P'JKTTOTl 

Uld)  = FU2(TT0T) 

IF<!pLAG>  40,40,59 

40  PR  I‘1T  41,L0*1  Cl)  ,L0NC2) 

41  F0OMflT(1x/i7x,"L0'lGTTUniNAL  TIME  RESPON  SE“/20X,  2A 1 0//5X  , "TI  h-  " , 7 X 
1"U~,8X,"4L°H  A",3X,"n",8X,"r-CTrA",7<,"3E,\9X,"0r"> 

GO  TO  60 

<50  PRINT  51,  L6T  Cl),  LIT  C2) 

SI  FOP.MAKix/ROX,  "LATERAL  Til!  RE  5P0N3£*V21< , 2A19/ /5X, "T HE",SX  , 


1**9£T  A~,7X  , “PV,  ox  , “R",  aX,"R-<I”,  3X  ,"0W“»9X» “OR") 

30  T3  30 

60  PRINT  SI, “TN IT  COMO",  CXC  I)  , I=!  , *)  , JIC 1)  , JIl  2) 

61  FOR’ AT  (1X,A9,?X,C9«5»  IX,  P9 . o ,?X,  =-R.  5 , IX , r 10. 5, 2 C2X,  Fq.  5 ) ) 
CALL  NATNULT  CAI,  AIAT,  AHj-ijS,*!) 

CALL  MATMULT  C4H,3-‘AT,H,M,'I,  ?) 

70  OKI)  = FUlCTTon 
UK?)  = FU2CTT3T) 

CALL  MATMULT CE4T,X,Gx,M,N,l) 

• CALL  HATN'JLT  CH,'JI,Hli,",2,l) 

CALL  MATA00CGX,H'J,X,M,1) 

TTOT=TTOT  *T 

Ic  Cl UNITS • EO  .0)  GO  TO  75 
00  75  1=1, « 

75.  XP<I)  = X (I)  *57.  295 
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Table  III  (cont.) 

U°<11  = UI(l >*57. 29* 

(JP<p)  = t)T  (?)»5’.  295 
IF  <TCLAG«  Lr.  1)  X°(1>=X<1) 

IFCFLflG.i.t.  3)  J=>  <?)-=•»!(  2) 

33  T0  30 

76  o:  ’ 7 1=1, M 

77  X=(I>  = X <I) 

U°  (1 ) = 'JI(l) 

U=(7)  = UI  (2) 

30  o?I‘iT  91,  TTpT,  (XD(I>,  1=1, *1)  , JP  (2) 

31  n IX ,F1  0.9, IX,  F9.5,1X,P9.  5,2X,F9.9,K,Fi0.5,2(2X,ir9.t;)> 
IFdPUP)  90,91,100 

C 

c*****»»  'IPOATS  LONGITUDINAL  PLOT  A99AY 

C 

90  IPLPT  =1 PL  0T+ 1 

IFCOLOT.  LT.  N=L3T)  GO  TO  35 

IPLOT  = o 

N°1  = N°1 +1 
TlCtPl)  = TTOT 
U(N31)  = XP<1) 

AL(',D1>  = *p(2) 

.Q(H»1>  = X°  (.3) 

TH  ('IP  1)  = XP(4) 

Oi<MPl>  =MO(l) 

9TCIP1)  = UP(  2) 

95  IP(TTOT.LT.riAX)  GO  TO  70 
GO  TO  110 
C 

C*******  "PQATc  LATERAL  PLOT  ApP.AY 
C 

100  I°LOT  = IPLOT+1 

IFCPLOT.LT.  JJPLOT)  r.O  TO  105 
IPLOT  = 0 
NP2  = NP2*1 

Tei'jp?)  = nor 

9i(NP2)  = X° (1) 

P(MPP)  = XP(2) 

<CNP2)  = X°  (3) 

PH  (NP2)  = XD(!») 

OH  ('IP?)  =UP(1) 

OPCIop)  =(JO(2) 

105  tc(TTOT.LT.TMAX)  GO  TO  70 
GO  TO  110 

110  IF(IFLAG)  131, 120 ,1 T0 
120  IrLt  G = 1 
GO  TO  10 


% 

M 

I 

.5 


C 

C*»**M*  °LOT  STATS  tf A^IAHLS S 

c 


130 

CALL 

OUTPLOT  (SP1.NP2) 

GO  TO 

1 

1M 

STOP 

"IHSUFC IOISNT  OATA  C49DS” 

150 

STO° 

~A  M1T9IX  IS  SINGULA?” 

160 

STOP 

.EHO 

129 
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Table  III  (cont.) 

SUO^OIJTIN  E fftTi^T  (AMA’VIAT  4 i"4T  , 1,  'I,  0 »T ) ^ 

0 1 Mr  MS  I OH  A*  AT  ( 1,  M)  , AINT  (!*»*•>»  A =01  ( ^ ♦ » *4,  4) 

01 M'  MCI  ON  TO  t4,4  ) , ASU  1(4 ,4 ) » "AT  ( M,  N)  , Ai  AT  ( N*N) 

imt"g:p  m,n, o 
00  ? 1 = 1*  -1 
00  *»  J=1,N 
I0(T,  J>=3 

AINTII,  J)  = A*<  AT  ( t t j)  »T 

00  T 1=1,  M 

I0<T,TJ=I0(I»T>*1 

00  " K=2,0 

OALL  lATp0H(A^AT,4P0H,H,M,><) 

00  4 1=1, - 

00  - J=l,'l 

G<I,  J)=(4°0H (T, J)»(T**K))/rS0T(<) 

CALL  1ATA3C<G,AIMT,AS'J«,1,N> 

00  • 1=1,  1 
00  " J=1,.H 
AINT(T,J)  <i  , j) 

00  1 1 = 1,  1 
00  5 J=l» N 
AEAT  (I,  J)  =AIMT(T,  J) 

EAT<T,J)  = Ai-AT(l,J)4.n(I,J) 

RETURN 

ENO 


InE  PATAOO 


7L/7L  03T=1 


= TN  4.5+-»l4 


SUBROUTINE  MATAn(A,O,0,i,'J) 

DIMENSION  A<M,N)  ,3<*Vi)  »C(N»N> 

00  1 1=1,1 

00  1 J=1,M 

1 C<I,J)=A(I,J)»0(I,J) 

RETJPN 

rsn 

‘INE  MATPOW  . 74/74  0PT=1 

SU^OIITIN"  *ATPOW(A,C,N»N,05 

DIMENSION  A<M,NJ ,3(4,4) ,C(M,N'  ,0(4,4) 

INTEGER  0 

IC0’INT=1 

00  1 1=1, 1 

00  1 J=  1 » N 

0(I,J)=A(I,J) 

1 3(1, JJ=A(I,J) 

3 00  ? 1=1,1 

00  ? J=1,N 

cex, j>=o 

00  ? K=l, N 

2 C<I,J)=C(T»J)*0(I,<)*0(<, J) 
IC'OJNT=IOOUMT*l 
IMICOUNT.GE.OJSO  TO  100 
00  V 1=1,1 

00  4 J=1,N 


FTN  4, 5*414 


oci,j»=c<i,jj 

50  TO  3 

CONTINUE 

RET-J^N 

E10 


f 


Table  III  (cont.) 


FUNCTION  Cflc 

int-str  m,« 

A = 1 

00  1 K = 1,N 
1 A = A*< 

FACT -A 
RIT')Pfi 
r.jn 

:n- 


r(*o 


7U.nL.  03T=1 


rTN  4.5+-.1-. 


S'JPPOUTIN*7  HATM'iUT(  A,  3 ,C  , M, M ,0) 

01  Mr  NSION  A ( H,  N)  , 3 ( “,  0)  , 0 ( H , 3)  , 0 (4,  4) 
IMTir.rP  M,N,  0 
00  1 1=1. H 
00  t J=1,H 

1 0<I,J)=A(I,J) 

00  ■>  1=1,  M 

O0  2 J=1,0 
0(1, J)=n 
00  2 K=1,N 

2 ctli  J)=C(I,J)+0(I,<)*9(K,J> 

RrT'JPfl 

E >10 

:NE  OIITP'.OT  74/74  0°T  =1 


PTH  4.5+41-. 


i 


suo3onrrsii  outplot  (ni,N2) 

COM‘<ON/L  NG/TinnO)  ,U  (305)  , 0 t?  0 5 > , AL  (305)  ,TH(305)  ,01(305)  ,0T  C*05) ; 
1 LON (2) 

COMMON/ UT  /T2(*Q3)  , p (305)  , ®.  (3 05) , 3£  (70 5)  , PH( 30  5)  , 3W  ( 30  = ) .O’ (305) 
1 UT(2) 

CAUL  PLOTS  (3  0) 

CALL  PLOT (1., -5.  ,-3) 

XF  Oil.  £0.9)  GO  TO  30 


C*******  LONG IT OOI NAL  V ARIA9LF S 
C 

CALL  SCAL  E (T 1, 4. , N1 , 1)  ' 

M=0 

C 

C *******  “LOT  OE,  OT 
C 

CALL  PLOT  (1,  , 1. , -3) 

10  CALL  AXIS  (0,  ,C.,  11HT  (3EC0N0  S)  , - U , 4,  , I) . , T 1 (Ml*  1)  ,U  ( Nl*2 ) ) 
CALL  AXIS(0. ,9.  , 1 2H0E  (023 R”3 ) , + 12 , 2.  , 90 . ,-29.  ,29.  > 

• CALL  AXIS (4.  ,0., 10HDT  (UNITS ) , -10 , 2. ,90., -10, ,10,) 

Ot(Ml+l,)  = -20. 

OE  <*!i+2)  = 29, 

OT(>!x+1)  = -10, 

OT(M1+2)  = 10, 

J1  = Nl'/io 

CALL  FLIN2 (T 1,02, -H t» l» Jl, 1) 

CALL  FLTNE (T 1,07  »-Nl, l, Jl» 2) 

CALL  SYM3  0L(  1,59,2.  ,,74, 129, 9.  ,-U 
CALL  SYM90L(1.90,2.,.07,  32,  0.,-l) 

CALL  SYM30L(  2.00,2.  ,.  21,  77,  0.  ,-l> 

CALL  SYP70L(  2. 21, 2.  ,.24,123,9.  ,-l) 

CALL  SYM30L(  2.42,2.  ,.  07,  47,0.,-l) 


IF(M.EQ.l)  GO  TO  20 

N«1 

CALL  PLOT  (G.  ,0.,-3) 
60  TO  10 
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Table  III  (cont. ) 

• C 

°L OT  flLP-U 
C 

20  CALL  DL  OT t-6.,3.0,-3) 

CALL  SCAL"(AL,2.  ,‘ll,l> 

CALL  AXTStO.  ,5.  ,!M  ,-L ,4.  , 0 . , T 1 (N1H  ) ,T  l (Ml*?) ) 

CALL  AXIS  tp.  ,<).,«*«( ''SORE'S*  , *3, 2.,  90.  ,AL(N1*1>,  AL(U1*2>  ) 
CAL'.  rL  TO"  (T  lj  At » •'!  i f t , 0,L) 

CALL  SYMJOLt?.,:’.  ,.23,  01,0. ,-l) 

C 


C*******  °L0T  TH£t4 

cj 

; CALL  PLOT  (*.  ,0,,-T) 

I call  SCAL£(T«,’.  ,*ll,l> 

CALL  AX  IS  ( P.  , 0.  , 1 H ,-1,4., 1.  .’•i(Nlfl)  ,Tim  + 2>) 

CALL  AXIS  (0,  ,0.,OH( O'SHEES)  , *-0, 2.,  0 0.  ,TH(‘U*1>  , TH(Nl+2>  ) 
CALL  PL  IN  E (T  1,TH,-M1,  1,  0 , L > 

CALL  SYM00L(,.,2*»«21»1 44,  0 , 1 > 


C 

C *******  nL OT  ') 

C 

IN;  OUT°LOT  74/74  0=>T=1  STN  4. 5*414 


CALL  t>LOT  C-S.,3.  0,-3) 

CALL  SCALK  U,2. ,Ni,l) 

CALL  AXIS (0,  ,0, , LH  ,- 1,4. , 0 . ,r 1 ( Nlf 1) , T1 (Ml*?) ) 

CALL  AXIS(0.  ,i.,90  (PT/sr~)  , ♦.?,?  -30.,  Om+1),  OtNl*2>) 

CALL  PLINKTl,  0,-Nl,l,  C,l> 

CALL  SYH30L(2.,2.,.21,  48,0. ,-l) 

C 

C*******  “LOT  o 
C 

CALL  PLOT  (P.  ,0.,-T) 

CALL  SCALK  0,2.  ,*lt,l) 

CALL  AXIS tO. ,0.,1H  ,T1(N1*1>  ,T1(N1*?>) 

CALL  AXIS ( 0.  ,0.,9U(OEG/3EC)  ,*0,2, ,90.,  0(01*1),  Q(*J1*2>) 

CALL  CLIN £ (T 1,  0,-Ml,t,  0,t) 

CALL  SYMOOLt 2. ,7.  ,. 21,  4i,0.,-1> 

CALL  SYH3  0L(-3.TS,2.3,,1S,2&HL0'I0ITUDIN4L  TIMS  RES=>0‘ISE , 0 . , ? 5) 
CALL  SYM90L(-2.5,2.  ISjLD'1  tl)  ,0.,10) 

CALL  SYMOOLt  -1, 0 , 2,  c-, . 15,L3N  (2)  ,0.,10) 

C 

c *******  lateral  variables 
C 

30  IF  ('ll.  EO,  0)  CALL  °L OT  ( 0.  ,5.  , -3) 

IF  04?, SO,  0)  G3  to  60 
CALL  SCALE(TP,4. ,N2,1> 

H*0 


ooo 


Table  III  (cont.) 


-*.  k"05/«’  v»’  > 


♦***••*  c-toT  op*  ok 


'*(1 


CM1.  PLOT  <’.  ,-S.  , -3) 

CMC  AXIS  <<>.  HMT  CS--CO*J"S)  , 

cat-  axis  (0.  ,1. ,q^(n"'PciS> , *0, 
0=  ri2  + l)  = -20  . 

<p  ei2*2)  = 2r. 

OW  (*1 3 ♦ 1 ) = “2n  . 

04(i’*2)=2"» 
j?  = tJ2/in 

CAUL  FLINc<T2*np,-N?,l»J2*l.) 


-11,4.,0.,T2(12*1>  , 

2. , 9n. ,-20. *20. ) 


t2C)2+2)  ) 


Oft LL  FLTN2<T2,0W,-N?,1,J2»?> 

CALC  SYHT0L(1.S9,2. 

CALL  SYMOK  t.  00 , 2.  *.  07,  *5,3.,-l) 
CALL  SYM30L(?.01, 0.  ,.?1,  7T»1.*_1> 
CALI  «YK10l(  ’•  "*1  »■’•*.  74,1’9  , n.  , ” 1> 
CALC  «Y^90L( **5 , 2.  ,.  32,  53,3. , -1> 
IFC-MTO.l)  GO  TO  30 
1=1 

CALL  °LCT  (P.  ,0.,-T) 

*0  "0  40 
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Table  III  (cont.) 

CALL  oL0T(-6.,3.1,-3) 

CAUL  SCALE  CP",  2.  ,'!2,D 

CALL  AXIS  (0.  ,3.  , 1H  ,-1,4, ,0.  ,T2(N2+1)  ,T20l2+2>) 

call  AXIS<3.  ,0.  , T-*(  fl"  <«rf  ?>  , +3, 2. ,30.  ,BE(M2+1)  ,-3S(N2  + 2>) 

CALL  FLlNE  02,33  ,-N?,  1,  ■',!) 

CALL  SYM30L(2. ,2.  ,.23,  32,0., -15 


Mi  OiJTPLOT 


74/74  0DT  =1 


•TN  4.5+414 


“LOT  PHI 


C 
C « 
C 


call  PLOT  <6.  ,3. , -3) 

CALL  SCAL £(P^,2.  ,M?,1) 

CALL  AXIS  (0.  1«  ,-1,4,, 0.  ,T20)2+1)  ,T2(M?+2>) 

CALL  AXIS  (0.  ,0.  , 9w<n£G=SE2)  , +a , 2. * 30.  »PHOJ2  + l)  , PH  ( *12  + 2)  > 
CALL  FLIN£(T'3,3H,-M2,  1,  0,1) 

CALL  SYM3  0L(2. ,2. ,.  23 , 140 , 0 . 1) 


'’LOT  o 


CALL 

CALL 

call 

CALL 


=LOT (-5. ,7. 0,-7) 

SCALE  ( 3,?.,N2,1) 

AXIS  ( 0,  ,0.  , 1H  ,*1,4,  ,0.  ,T ?((I2+1)  ,T2  012+2)) 

AXIS  (0.  ,0.  ,94<M£S/srC)  , +3, 2. ,30.  , °0)2+l),  °<N2+2)) 
CALL  FLTN * (T  2,  3,-H2,l,  0,1) 

CALL  SYH30H  7., 2.  ,.21,  4?,0.,-l> 


C 

c 


3L1T 


60 

C 


CALL  PLOT  (6.  ,0., -3) 

CALL  SCAL£<  9,2.,»I2,1> 

CALL  AXIS  (0. ,T.,1H  ,-1,4. ,0. ,T2(N2  + 1) ,T2(M2  + 2) ) 

CALL  AXIS(0. ,0.,9H(0£S/S-C) ,+P,2.,30.,  9(M2+1),  9CN2+2)) 
CALL  CLIN  £ (T  2,  9,-M7,t,  0,1) 

CALL  SYH9  0H2.  ,7.,»21,  45,0. ,-l) 

CALL  SYM3OL(-2.30,2.8,.13,2HL5r£9SL  TIME  RESPONSE, 0 . , 2 1 > 
CALL  SYM30L( -2.3 , 2. 5, . 15,1  AT  (!) , 0.,  10) 

CALL  SYM30L(-1.0,’.5,.1S,LATC7)  ,0.,10) 

CALL  PLOTS <1  DUMMY) 


RETURN 


EMO 


ON  Fill 


74/74  OPT  =1 


FTN  4.5+414 


ON  FU2 


FUNCTION  Flll-(TTOT) 

EL EVA TOP/ AIL  CRON  zero  inpjt 
FU1  = 0. 

RETURN 

ENO 

74/74  0®T=1 


FTN  4.5+414 


FUNCTION  F!J2  (TTCT) 

THRCTTLF/ RUT  3SR  IM=ULSC  IM=UT 
COMMON/UFUNC /T, I FLAG 
TIMP  = 0. 

FU2  » 0. 

IFIIFLAG.LE.  O.A'JO.TTOT.GE.  TIMP  ) 
IF(tFLAG.  r,+,  n.  AMO.TTOT.GE.  TIMPP) 
IF <TTOT.G£,T IM=  + T ) FU7  = 0. 

RFTII9M 


FJ2 

FJ2 


= 5. 

= .175 
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Appendix  D 

FLEXSTAB  Evaluation,  Flow  Charts  and  Input  Listings 


This  Appendix  evaluates  the  FLEXSTAB  programs  used  in 
this  report  by  discussing  their  overall  effectiveness  in 
meeting  the  requirements  of  the  previous  sections.  Flow 
charts  and  input  listings  are  also  provided. 

Generally,  the  use  of  the  programs  needed  in  this  ana- 
lysis is  simplified  by  a wexl  written  user's  manual.  Many 
of  the  programs  can  be  executed  successfully  only  after 
one  or  two  trial  and  error  attempts. 

The  Geometry  Definition  (GD)  program  of  FLEXSTAB  is 
the  easiest  to  understand  and  to  use.  However,  some  minor 
improvements  in  the  plotting  of  the  model  would  be  very 
helpful.  For  instance,  the  fuselage  camber  definition 
needs  to  be  plotted  in  order  to  give  the  user  a better 
physical  picture  of  his  results.  Also,  as  a minor  change, 
the  slender  body  fuselage  definition  plot  should  be  shown 
closed  at  both  ends.  There  is  some  possibility,  in  pre- 
paring the  data  for  the  fuselage  definition,  in  which  the 
user  may  in  error  leave  the  nose  and  tail  open.  If  not 
desired,  this  creates  erroneous  calculations  of  the  model's 
stability  and  control  characteristics.  These  problems  are 
avoided  by  using  the  program  listing,  but  a more  detailed 
plot  of  the  airplane  is  by  far  more  efficient  in  spotting 
errors . 
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The  Aerodynamic  Influence  Coefficient  (AIC)  program 
of  FLEXSTAB  is,  internally,  the  most  complicated.  This 
program  uses  the  GD  data  to  define  the  model  aerodynamically . 
However,  this  is  done  for  only  one  Mach  number  each  time  the 
program  is  executed.  Due  to  the  high  time  and  cost  of  exe- 
cuting AIC , convenience  in  changing  the  Mach  number  for 
the  model  is  not  available.  This  is  a possible  improvement 
area  needed  in  this  program. 

The  Internal  Structural  Influence  Coefficient  (ISIC) 
program  is  far  more  complex  than  the  GD  program  in 
defining  the  physical  characteristics  of  the  model.  It  is 
also  the  most  difficult  to  use,  particularly  in  trouble- 
shooting errors  in  input  data.  For  example,  even  though  all 
the  structural  information  input  into  ISIC  must  be  defined 
for  all  slender  and  thin  bodies  defining  the  model,  the 
program  will  execute  successfully  even  when  this  requirement 
is  not  satisfied.  Needless  to  say,  this  could  have  serious 
effects  on  subsequent  programs  where  the  results  of  ISIC 
are  used. 

The  Stability  Derivatives  and  Static  Stability  (SDSS) 
program  uses  the  results  of  both  GD  and  ISIC  in  its  calcu- 
lations. This  program  does  not  create  significant  problems 
in  preparing  input  data  but  soi  - improvement  is  needed  in 
clarifying  the  program’s  output. 

The  static  and  dynamic  stability  derivatives  that 
are  listed  by  SDSS  are  in  the  stability  axes  system  of 
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the  model.  However,  no  where  on  the  listing  does  it 
indicate  this.  This  creates  problems  if  the  data  is  sub- 
sequently used  for  a dynamic  analysis,  as  it  was  in  this 
report,  and  the  incorrect  body  axes  system  is  used.  How- 
ever, this  is  only  a problem  when  a dynamic  analysis  uses 
programs  external  to  FLEXSTAB. 

The  stability  derivatives  on  the  SDSS  listing  are 
generally  quite  clear  and  are  shown  in  a standard  format. 

The  one  exception,  however,  is  in  the  speed  stability  deri- 
vatives. As  listed  by  SDSS,  CL(U),  CD(U)  and  CM(U)  are 
not  the  same  as  Clu,  Cdu  and  Cma.  The  FLEXSTAB  values  are: 

CL (U)  = Cl£  + 2Cli 

CD^U)  = CD~  + 2CDi 

CM (U)  = Cmu  + 2Cmi 

This  is  not  clear  on  the  output  listing  and  some  improvement 
here  is  also  needed. 

in  tabulating  the  rigid  and  elastic  model  FLEXSTAB 
data  with  flight  test  values  in  section  II,  there  is  good 
correlation  seen  in  those  values  compared.  However,  refer- 
ring to  Appendix  E,  there  are  a number  of  other  FLEXSTAB 
derivatives  that  are  shown  but  not  used  in  this  report. 
Generally,  these  values  are  unsteady  derivatives  and  are 
not  easy  ho  interpret .physically.  They  are  assumed  to 
have  little  or  no  impact  on  the  research  objectives  and 
results  of  this  study.  Also  lacking  flight  test  or  wind 
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tunnel  data  for  making  comparisons  precludes  any  attempts 
to  validate  them.  Note  that  the  derivatives  used  in  this 
report  are  those  included  in  the  aircraft  equations  of 
motion  discu.ssed  in  Appendix  A. 

FLEXSTAB,  with  only  a few  exceptions,  is  seen  to  be 
an  accurate  and  effective  analytical  tool.  It  appears  to 
do  well  in  predicting  rigid  and  elastic  aircraft  character- 
istics withir  .e  limitations  of  nonviscous  linear  potential 
flow  theory. 

Figures  1 and  2 of  this  Appendix  display  the  two 
general  program  sections  of  FLEXSTAB.  The  cross  hatched 
sections  shown  were  not  used.  Sample  input  listings  for 
GD,  ISIC  and  SDSS  are  shown  in  Tables  I,  II  and  III,  respec- 
tively. 
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FIGURE  1 -FLEXSTAB  2.01.00:  AIRPLANE  DEFINITION  SECTION 
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Table  II  (cont'd) 


! 

! 

1 

1 

o 

1 

1 

• a 1 

O 

o • : 

o 

• 

1 

to  o i 

*4 

i 

i 

\D 

in  * 

*4  fO 

k 

in 

«*  M , 

N • 

’ 1 

eg  iD 

• «r4 

l 

• • 

*4  in 

«o  in 

to 

vo  in 

40  *4 

*4  PO  PO  PO 

H H ^ 

Ul  Ui  UJ  UJ 

sO 

O O C 

> o 

in  j-. 

*4  to 

44  o c 

> a 

j- 

in  in 

• • • • 

40  O 

*4  • 

fv.  <H  -«H 

o*  in 

• 0 

• • 

*4  N. 

1 { 

PO  40  . 

eg  eg 

•*4  cr 

• *4 

*4  PO  PO  PO 

« *4 

if) 

UJ  UJ  UJ  UJ 

in 

sO 

o o c 

> o 

▼4  K 

*4  fO 

«0  o c 

1 o 

n*  eg 

in  cr 

• • • • 

-*•  in 

«o  • 

«D  wi-r 

• *< 

O O'  , 

cn  • eg 

• • ■ 

o*-na 

in  k i 

• Z O'  eg 

1 

40  PO  » 

PO  lit  K-  *4 

eg  fo  to  to 

X 

♦4Hf 

1 *4 

o 

UJ  UJ  u 

1 iu 

UJ  >c 

O D C 

1 o 

W4« 

•*•4  0 0 0 

▼4 

in  o 

uiuiuujuitijujajuiui^juj 

ODOOOODn  OOOUN 

oinegvO^njtoroo'cotT'Po 


^ M O'  O'  N -3*  -T  fOWr^H 


N 1 

eg 

i • 

eg  s 

1 

«D  O' 

- i i 

eg  *4 

• • 

; ‘ 

m to 

« i 

vO  *4 

in 

»4  H 
*4  ^4 

UJ  UJ  - 
o in 

N.  O' 

K to 

vO  40 

• • 

in  n. 

o o 

• • 

o to 

> 

*4  *4 

i 

UiUJldUJUJIiJUJUiUiiJUJ'd 

o^roo^roNiOoaorao 

ou)j>fOK'^®(owiroco 


UiluUJUiUJUJUJUiUJuJ'jJU 

CJQO'^vCHOOUOOQ 

oU'in^vDN^^SWOJN 


i 

Ul 

• • • 

o 

n-  ^ <t 

CO 

eg  o'  ^ 

• 

1/1  H 4 O 

CO 

J 

cn 

•»  J4 

eg  uj  *4  , 

UJ 

< 

p-  cn  ! 

o 

X 

4i  cn  , 

4t  ' 

1 

1 

Z \D  UJ  { 

X 

1 

o 

in  in  eg 

M4M2 

M 

UJ 

\D  eg 

o in  o u 

a 

O'  in  <c 

a:  in  • u. 

o 

X O 

eg  in  ^4  o' 

o 

o • n»  *h 

a 

0 0 0 0 3 • 

to  • m * 1 

• 

OvCsDH-OOOoOo 

«i 

• • 

• — j -3" 

• o o'  a 

eg 

u 

T 

i 

0 

0 

0 

a 

2 

3 

6 

$ 

1 

$ 

6 

1 

$ 

0 

0 

0 

0 

0 

0 

u>  to 

cn 

o o . 

• • 

O P- 

vO  N 

• z 

<0  N. 

^ UJ 

X 

*4  *4 

o 

*4  *4  . 

•* 

UJ 

UJ  UJ 

eg 

in 

o o ; 

»c  in 

o m 

o in 

u.  < 

• • 

cn 

• eg 

O f- 

*4  a <n 

UJ 

O'  O' 

O 41 

UJ 

cn 

N • 

• cn  o 

cn 

cn 

*4  pw 

PO  UJ 

CO 

4t 

p-  cn 

<r 

X 

1 < cn 

X 

eg 

Z UJ 

o 

■H  PO 

in  in  z 

o 

UJ 

iD 

a <©  u. 

UJ 

Q. 

C"Oh 

a • u. 

a 

X 

o • *4  or 

O J H 

X 

• o eg  «x 

0 

0 

c 

T 

0 

0 

u 

0 

3 

«4  *4  • OL 

• cj  «r  in 

• • — i • 

W4(\J  vfl  WfU  tftr(  WOOM  4 

> .3  in  tX> 


N®fl*OT4W/‘Jl^U)N«(rOv4WW«tU'^K«HJ»o^MW^U\iON. 

-j-J’-j-inininininininininin^gDvouDvDvOvOovc^r-N.r'-r^N-N.n.p^ 


c <J>  o 
N K ® 


t< a g ■ex^+H-.r&irzr At .-; -•<•  W u ' 


sfecMS 


IC“135A-NO~  WINGLETS - - - - - PAGE 


• • • 

• cp  cm 

K *4  & 
«o  *4  cm  ro 


UN  • 

• • \D 

• in  4 .4 

MOK  «4 


• MO  J 
4*  CNJ  « O 
n0  *4  <H  ro 


UN  in  ♦ 

• • fO 

• K-  4 «H 


C\l<7'O'O'(7'C'<TO'0r'C'<J'0r' 

*400000000000 

UiUlliJUiUiUiUJUiUiLdaiUj 

OQ  0^03  ITiOol^aOQ 
• « • • •ir^WHBOO 

•4  M O « W • • ••  • • • 

cvjcvj*4w<rr^u>inro(Nj^ 


O Ul 

■ 

V* 

Ul  O <VI  UN  *4 

X 

*4  * 

a 

X 

NO'O'tro'O'cro'u'O'O'cr 

o 

H 

Ul 

o 

*400000000000 

Ul 

Ul  1 

> 

UJ 

UiUJUiUJUJUJUtUUitJUjUJ 

CO 

O 1 

in 

OQO^  MSOU'Otnoa 

UN  1 

i 

U_ 

4*  • 

«t 

ul 

< *4  04  O s.  ro  O 

1 

o 

►-  O tO 

l- 

o 

• « • >-  CVJOJ*4*4*4<OvDinrOf\i*4</5 

c 

< Ul 

< o 

• 4 4 ro  < 

IU 

» 10 

a in  . 

c • 

1/5  UN  O U>  05  a 

in 

. J 

r*  uj 

in 

o 

m 

*- 

in  < 

V 

H* 

^ : 

< 

in  x 

a 

in 

X 

2 3UJ 

o 

z 

ui  ! 

i 

4* 

w in  z o 

CO 

M 

Z i 

o 

m 

Ot4  1l  UJ 

Q 

u ; 

Ul 

NO  *- 

o' 

• u.  a. 

z 

Of 

u 

a 

4-  C3Bf 

OvOM  31 

O'  M 

0 1 

O 

fO  IH 

X 

• CoO'Ot-oDo 

ro  X o 

< 

• o 

• • • r-  000000000003 

O £X  • O 4 (/>  • _J  • • ►— 

C\J  V*  I 


i • o in  fo 

» O W 3 O' 


»—  000000000003  *000^00 
WOOOOOOOOOOOOH^OOtyHOW 


*4njK>4inNCn-«o(7*o  <4(\jf04inv0^«oa'o*4<\jw,*-in\o^«>cr'« 
co^«octf(O0«o(TO'(r(ro>()'(r(ra>o'ooaooooooo< 


Hi.w^intAKco'o 

HT<v4r4r4T4<HHHM 


Table  III 
SDSS 


\ 


</i  c/>  c 

r 

a 

aJ  © 

c 

o 

S. 

© _i 

<r 

• 

1.0 

♦-*  J «- 

-j  r 

a 

• » » 

►-»*.© 

UI 

* 

<x  o o 

a 

V)  » 

w k 

o 

a * 

wzu. 

X 

or  • 

«r 

a 

«x  * 

o o a 

c a <t 

o * 

z z 

*-i  •«.  >— 

in 

» 

< • 

O’liU 

fw 

o 

• 

«*  » 

c c 

►4  CL  *-* 

• 

►-  «► 

to  z Z 

or  © <x 

• 

o 

< * 

uj  *1  UJ 

© • 

> \ 

»- 

t 

w 

to 

fr-  • 

K W ♦-» 

•% 

vD 

3 • 

«I  UIl 

5: 

© 

o 

a *► 

> -J  «r 

x 

« 

• 

z * 

w o 

z ♦-* 

© 

o o 

cv 

H * 

Q Z * tr 

o 

C' 

o \C 

• • 

x 

» 

.1  H I Z 

M 

• 

• if' 

IP  ri 

j- 

u.  «*■ 

CX<Ck J 

o 

CD  «4 

~4  » 

44 

o % 

x *-*  o 

» 

t-  t 

»-  »-  z «t 

t">  »l 

s 

<r 

to  % 

hhi.  C 

IJ  t C 

o 

h- 

J JtiHcr 

*•4  u u z v 

QC  < U4  M o 
<r  • c:  a r* 
H K .t  «*’  M T • 


• C C.  OL 

vD  U o Ui 

C*  -J  «4  it  4000*JC-00©0C.»C:0  coo 

* UN  C*  -I  X n>  ••«•«•  • * ••  • « 3 •• 


V) 

t)  *4  a 

U.  C\J 

<2 

•— » 

O 

CVJ 

t- 

44 

*-i 

O 

< f 44 

v4 

41 

ri 

^4 

4-» 

*4 

4-4 

44 

<r 

44 

44 

• j 

1 ►-  to 

IJ 

< 

►- 

z 

-J 

O 

«uz 

rv 

C 

Z 

c 

C X u -1 

a 

Hi 

o 

C -J 

z 

u.  y < 

c 

o 

• 

t- 

c 

CJ 

X !-• 

z 

•>* 

►4 

13  o 

Ul 

•J 

• 

to 

«3 

1-  < 

O 

a' 

w r • •%: 

.J 

a 

t- 

o 

if 

z 

X 

t-l  h* 

O' 

U‘ 

«->  c 

e 

10  5 « 1-  z 

u. 

<x 

It 

CD 

m 

(** 

u; 

►4 

© 

O 

O 

O Ul 

o 

QC 

-4  ua  u. 

© 

>- 

V 

m 

u 

<rr 

j: 

f • 

in 

a* 

ZOO: 

a t3 

o 

O 

o 

O 

O 

O 

u 

• 

• 

4 -J 

o 

z 

• © 

i. 

►- 

< 

C C «l  O 

o 

c 

z 

- 

• 

if 

i— 

• 

X 

% 

t- 

c • c 

• • 

• 

• 

• 

• 

• 

• 

• 

' * 

Of 

4 

*-• 

if  z 

• 

*4 

z 

1 

O 

*»  C O V*  O 

o 

►4 

%M 

=3 

r*"' 

♦h 

L_ 

44  M 

-J  **  X 

»4  (V, 

fo 

Jf 

m 

vC 

n- 

X 

cr 

•4 

44  *t 

H 

CVi  O' 

44 

«e 

r* 

• 

o 

1 

2 

3 

4 

Li'tCN 

X <r 

O 

4“  If. 

ifl 

K 

« «r  o 

**  Pi 

K* 

-*  in 

K 

or 

© 

t J 

•-4 

0;  O' 

^ ir 

cO  K 

O 

• 

z 

r< 

•4 

4» 

t4 

♦4 

44 

n 

tI 

•4  r*  *J  \ rj  f J 

rj 

CM  CM 

CM  CM 

rj 

f'' 

ro 

1-0 

fO 

fc,  fO 

po 

a 

it' 

•»l  C3  O 

I . - 


i tfo 
• * .«  • 

44  > K* 


VI 


153 


a 


Table  I 

FLEXSTAB  Rigid  Model  Data  (cont'd) 
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FLEXSTAB  Rigid  Model  Data  (cont'd 
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Tabic  II  ■ 

FLEXSTAB  Rigid  Short  Tail  Model  Data  (Flight  Condition  2A) 

*Short  Tail* 

No  Winglcts  With  Winglcts  With  Winglcts 

(Rigid) (Rigid) (Rigid) 


Mach 

.77 

.77 

.77 

h (ft) 

28,500 

28,500 

28,500 

V(fps) 

771.47 

771.47 

771.47 

q(psf) 

279.98 

279.98 

279,98 

Y (deg) 

0.0 

0.0 

0.0 

0B  (deg) 

-.89 

-.941 

-.941 

“w  (deg) 

1.11 

1.059 

1,059 

e.g.c 

.321 

.321 

.321 

- 

:XX  * 106 

(sl-ft2) 

2.05 

2.05 

2.05 

V x 106 

(sl-ft2) 

2.46 

2.46 

2.46 

I x 106 
zz 

(sl-ft2) 

4.36 

4.36 

4.36 

I x 106 
xz 

(sl-ft2) 

- 

- 

- 

Weight  (lbs) 

130000.0 

130000,0 

130000.0 

Tl  (lbs) 

2600.0 

2360.0 

2360.0 

6strim  (deg-) 

-.83 

-1,26 

-1.26 

♦Mass  moments  of  inertia  wrt  principal  axis  system, 
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Table  II 


FLKXSTAB  Riiiid 

Short  Tail  Model 

Data  (Flight  Condition 

, 2A)  - coat'd 

I 

I 

Lateral 

1 

*Short  Tail* 

1 

No  Winglets 

With  Winglets 

With  Winglets 

ii 

v 

t 

(Rigid) 

(Rigid) 

(Rigid) 

% 

i 

cYg(i/deg) 

-.0116 

-.0134 

-.0111 

V 

5 

| 

% 

c,6  (1/deg) 

-.0026 

-.0034 

-.0029 

1 

1 

CN3(1/deg) 

.0032 

,0034 

,0024 

I 

CYp(l/rad) 

-.1786 

-.3218 

-.2976 

C^pCl/rad) 

-.4959 

-.5577 

-.5533 

s 

V1/rad) 

.0135 

,0344 

.0233 

1 

Cyr(l/rad) 

.6615 

,7532 

,6192 

1 

CZr (1/rad) 

.1892 

.2222 

.1975 

I 

CNrd/rad) 

-.2677 

-.2851 

-.2237 

1 

I 

1 

CY6A(1/deg) 

1.31  x 10~4 

,94  x 10-4 

,34  x 10'4 

c,6Aa/deg) 

-5.74  x 10~4 

-5,9  x 10-4 

-6.01  x 10-4 

CN6A^1/,deg') 

-1.32  x 10“4 

-1,27  x 10"4 

-.99  x 10"4 

'% 

CY6Rd/deg) 

.0076 

,0075 

,0066 

1 

CS,3R(1/deg) 

.00084 

,00084 

,00066 

| 

1 

cN6Rd/deg) 

-.0034 

-.0034 

-.0030 

“5 

2ft 

1 

jy 

Cygd/rad) 

-.1784 

.1746 

-.1988 

3 

1 

(1/rad) 

-.0022 

,0030 

.0016 

V 

a 

% 

CNgd/rad) 

.0036 

.0036 

.0104 

h 

GypCl/ra*^) 

-.155 

-.133 

-.1108 

, j 

C^-U/i-ad) 

-.0097 

.0055 

,0100 

CNp(1/ra<i) 

.0669 

.0637 

.0532 

1 

•i! 


Table  II 

FLDXSTAB  R ! "id  Short  Tail  Model  Data  (Flight  Condition  2A)  - eont'd 


No  Winglets 
(Rigid) 


With  Winglets 
(Rigid) 


♦Short  Tail* 
With  Winglets 
(Rigid) 


CYi>(l/rad) 

C^(l/rad) 

CN-(l/rad) 


-.0171 

.0022 

-.0192 


-.0185 

-.0001 

-.0193 


.0028 

,0019 

-.0275 


Short  Period 


Phugoid 


Dynamic  Modes  Comparison 
♦Longitudinal* 


Co  (rad/sec) 
n 

2.37 

2.65 

2.65 

Period  (sec) 

3,67 

3.03 

3.03 

,691 

.621 

.621 

^s.p. 

T%  (sec) 

.423 

.421 

.421 

un (rad/sec) 

.061 

.060 

.060 

Period  (sec) 

103.74 

104.53 

104.53 

?P 

.010 

.011 

.011 

T^  (sec) 

1135.44 

1044,84 

1044.84 

♦Lateral* 

Dutch  Roll 

o)^  (rad/sec) 

1.93 

1,98 

1.66 

Period  (sec) 

3.3 

3.22 

3 , 83 

?DR 

.153 

.168 

.158 

(sec) 

2.36 

2.09 

2.65 

Tabic  H 


FLHXSTAB  R i i d Short  Tail  Model  f)afa  (FI  i 

No  Wing lets  With 

(Rigid) (Rigi 


675.68  470,4 

203.4  141,6 

Roll 

TVlo  (sec)  1.242  1.113 

T \ (sec)  ,374  .335 


Spiral 

T1/ iQ  (sec) 
T\  (sec) 
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gi»t  Condi  t ion  2A)  - cont'd 


Winglets 

d) 


*Short  1’ail* 
With  Winglets 
(Rigid) 
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Table  III 

FLEXSTAB  HI— tic  Model  Data  (cont'd) 


Static  Stability  Parameters 


CRUISE  1 

No  Winglets 

With  Winglets 

Static  Margin 

(Rigid) 

.2235 

.2572 

(Elastic) 

.1709 

.1830 

Neutral  Point 

(Rigid) 

.4534 

.4893 

(Elastic) 

.4008 

,4152 

Neutral  Point  Shift 

-.0525 

-.0741 

(Elastic-Rigid) 

Maneuver  Point 

(Rigid) 

.4779 

.5145 

(Elastic) 

.4226 

,4372 

Maneuver  Margin 

(Rigid) 

.2480 

.2923 

(Elastic) 

.1927 

.2051 

Longitudinal  Control 
Per  G for  Steady  Pull-Up 
(Deg/unit  load  factor) 

(Rigid) 

-3.4111 

-3.9361 

(Elastic) 

'2.8686 

-3.0590 

Stick  Speed  Stability 

.0091 

,0100 

(Deg/ft/sec) 
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Table  I IT 

p rnxSTMi  Elastic  Model  Data  (cont.  O 


Static  Stability  Parameters, 


CRUISE 


No  Kinglets 


With  Wing lets 


Static  Margin 
Neutral  Point 


(Rigid) 

(Elastic) 


(Rigid) 

(Elastic) 


Neutral  Point  Shift 
(Elastic-Rigid) 


Maneuver  Point 


(Rigid) 

(Elastic) 


Maneuver  Margin  (Rigid) 

(Elastic) 

Longitudinal  Control 
Per  G for  Steady  Pull-up 
(Deg/unit  load  factor) 

(Rigid) 

(Elastic) 

Stick  Speed  Stability 
(Deg/ft/sec) 


.2318 

.2052 

.4617 

.4351 

-.0266 


.4743 

.4470 

,2444 

.2171 


.2684 

.2290 

,5005 

,4611 

-.0394 


,5135 

,4733 

,2814 

,2411 


-6,9187 

-6.3949 

.0205 


-8,0827 

-7,1555 

,0229 
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FLEXSTAB 

Static 

Tabic  III 

Elastic  Model  Data  (cont'd) 
Stability  Parameters 

CRUISE  3 

No  Winglets 

With  Winglets 

*? 

'% 

Static  Margin 

(Rigid) 

,2211 

,2507 

?; 

s* 

L.V 

(Elastic) 

,1870 

,2030 

1 

if 

Neutral  Point 

(Rigid) 

.4510 

.4828 

g* 

✓ 

(Elastic) 

,4169 

.4361 

sr 

1 

Si  . 

Neutral  Point  Shift 

-,0341 

r-,0467 

| ' 
§S* 

(Elastic-Rigid) 

Sr 

| 

Maneuver  Point 

(Rigid) 

.4928 

.5257 

? 

fc 

(Elastic) 

.4558 

,4756 

s' 

| 

fe 

Maneuver  Margin 

(Rigid) 

.2629 

,2936 

1 

(Elastic) 

.2259 

,2435 

1 

Longitudinal  Control 

| 

Per  G for  Steady  Pull-up 

§■ 

(Deg/unit/ load  factor) 

i j 

(Rigid) 

-5,9830 

-6,7604 

1 : 

(.Elastic) 

-5,3997 

-5.8477 

fe 

i 

Stick  Speed  Stability 

.0214 

,0234 

S 

1 

$ 

(Deg/ft/sec) 
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Table  III 

FLEXSTAB  Elastic  Model  Data  (cont’d) 


Table  III 

FLEXSTAB  Elastic  'todel  Data  (cont'd) 
Lateral  Data  CRUISE  3 


No  Winglets 
Rigid 

Elastic 

With  Winglets 
Rigid 

Elastic 

Cyg  (1/deg) 

-.0105 

-.0097 

-.0120 

-.0110 

C£B  (1/de^) 

-.0020 

-.0018 

-.0025 

-.0020 

CN6  (1/deg) 

.0030 

.0026 

.0032 

,0028 

CYp  (1/rad) 

-.1626 

-.1537 

,2790 

-.2433 

Cep  a/rad) 

-.4364 

-.3612 

-.4833 

-.3803 

V1/rad) 

.0277 

.0221 

.0441 

.0316 

CYr  (1/rad) 

,6414 

.5952 

.7319 

.6717 

C£r  (1/rad) 

,2264 

,1931 

.2530 

.2004 

CNr  C1/rad) 

-.2798 

-.2586 

-.2981 

- . 2730 

CY6A  (1/deg) 

.000095 

, 000068 

,000071 

.000059 

cm  '1/de« 

-.000486 

-.000412 

-.000495 

-.000416 

CNSA  (»/<*•*) 

-.000080 

-.000070 

-.000077 

-.000070 

SfiR  (1/deg) 

.0064 

,0052 

.0064 

,0052 

c»«  Cl/deg) 

,0003 

,0003 

.0003 

,0003 

CN6R 

-.0030 

-.0025 

-.0030 

- ,0025 

Cy*  (1/rad) 

-.1748 

-.1678 

-.1735 

-.1666 

C££  (1/rad) 

-.0010 

-.0001 

.0015 

,0014 

CNj  (1/rad) 

.0033 

.0001 

,0032 

,0000 

(1/rad) 

-.1200 

-.1129 

-.1057 

-,0°93 

ctp  «/rad) 

-.0085 

-.0068 

-.0001 

- . 0o^ 

% CVxad) 

.0544 

,0511 

,0512 

.0482 
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Table  III 

FLEXSTAB  Elastic  Model  Data  (cont^cQ 
Dynamic  Modes  Comparison 
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